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ABSTRACT 


This  thesis  is  concerned  with  several  aspects  of  the  role  played 
by  spin  in  general  relativity,  more  specifically,  in  particle  equations  of 
motion  and  in  field  equations  for  spinning  matter.  These  equations  are 
derived  from  variational  principles,  where  the  Lagrangian  has  some  depen¬ 
dence  on  terms  associated  with  spin,  in  addition  to  the  usual  dependence 
on  the  metric,  external  fields,  etc. 

Chapter  One  presents  a  general  discussion  of  the  thesis  and 
introduces  the  notation  used.  A  more  detailed  introduction  to  notation 
appears  in  Appendices  One  and  Two,  which  discuss  invariance  identities  sat¬ 
isfied  by  a  scalar  Lagrangian,  these  having  important  use  in  Chapters  Two 
and  Three. 

In  Chapter  Two,  covariant  equations  of  motion  for  a  single  spin¬ 
ning  particle  are  derived,  and  the  special  case  of  motion  in  an  external 
electromagnetic  field  is  discussed. 

Einstein's  gravitational  equations  for  a  medium  with  intrinsic 
spin  are  examined  in  Chapter  Three.  The  form  of  the  total  energy  momentum 
tensor  for  such  a  medium  is  found  and  it  is  shown  how  the  equations  of 
Chapter  Two  may  also  be  recovered  from  the  four  dimensional  action  integral. 

Chapter  Four  generalises  the  work  of  Chapter  Three,  an  action 
integral  depending  on  any  number  of  derivatives  of  the  fields  is  considered, 
and  finally,  the  derivation  of  particle  equations  of  motion  from  the  Dirac 
equation  is  discussed  in  Chapter  Five. 
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CHAPTER  ONE 


GENERAL  INTRODUCTION  AND  NOTATION 


The  main  concern  of  this  thesis  is  an  examination,  in  the 
general  relativistic  framework,  of  the  field  laws  governing  fields  gen¬ 
erated  by  matter,  where  the  matter  is  assumed  to  possess  intrinsic  spin. 

By  field  laws  we  mean  the  gravitational  field  equations  and  equations  for 
any  other  fields  present,  such  as  the  electromagnetic  field. 

For  a  medium  with  no  internal  spin  the  Einstein  equations 
G^V  =  T^V  are  the  usually  adopted  gravitational  equations,  coupling  the 
Einstein  tensor  to  the  total,  symmetric,  conserved  (i.e.  T^V | ^  =  0  )  , 
energy  momentum  tensor  T^  of  the  medium.  In  the  spinning  case  one  may 
either  keep  the  Einstein  equations  and  examine  the  constituent  parts  of 
the  energy  momentum  tensor,  to  find  what  terms  attributable  to  material 
spin  appear  in  T^  ,  or  one  can  try  to  modify  or  generalize  the  Einstein 
equations  themselves. 

The  latter  approach  has  led  to  several  new  theories  [1].  There 
is  the  generalisation  of  Einstein's  theory,  originally  due  to  Elie  Cartan, 
which  takes  space  time  to  be  a  four  dimensional  differentiable  manifold 
with  a  metric  tensor  and  affine  connection  compatible  with  the  metric,  but 
with  the  connection  no  longer  taken  to  be  symmetric.  The  antisymmetric 
part  of  the  connection  is  then  related  to  the  density  of  intrinsic  spin. 
Another  idea  is  the  suggestion  of  Sciama  [2],  that  the  gravitational  field 
produced  by  particles  with  spin  should  be  described  by  a  non-symmetric 
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2. 


metric  tensor  g^V  ,  the  energy  momentum  tensor  thus,  in  general,  being 
also  non-symmetric . 


In  contrast  to  the  above,  the  approach  taken  in  this  thesis  is 
to  retain  the  orthodox  theory  of  a  Riemannian  space  time  with  a  symmetric 
metric  tensor  and  symmetric  connection.  We  find  the  form  of  the  total, 
symmetric,  conserved,  energy  momentum  tensor,  a  form  containing  terms  due 
to  intrinsic  spin  and  one  that  reduces  to  the  usual  energy  momentum  ten¬ 
sor  in  the  absence  of  spin.  Variational  principles  are  adopted,  the 
contribution  from  internal  spin  stemming  from  the  dependence  of  the  Lagran- 
gian  upon  a  tetrad  field  ev  ' (x)  describing  particle  spin. 

Using  a  statistical  approach,  Werner  Israel  [3]  has  shown  that 
the  gravitational  equations  should  take  the  form 


-(8tt)  1G°^ 


=  Ja0  =  T°^ 


+  \  (saY3 


+  s 


$ya 


-  sa3Y) 


\y 


(i.i) 


for  a  medium  with  internal  spin  acting  as  source  of  electromagnetic  and 

cx  y  6  cx  3 

gravitational  fields.  Here  S  is  the  material  spin  flux  and  T  is 

the  energy  momentum  tensor  for  non-spinning  media  (the  sum  of  a  material 

cy  R 

and  an  electromagnetic  tensor).  T  is  neither  symmetric  nor  conserved, 

R 

only  the  total  energy  momentum  tensor  J  has  these  properties.  It  will 
be  seen  that  the  Lagrangian  approach  gives  exactly  the  same  form  for  the 
total  tensor  J  as  found  in  [3]. 

However,  the  first  part  of  the  thesis  is  concerned  with  the 
derivation,  from  an  action  principle,  of  equations  of  motion  for  a  spinning 
particle  moving  in  given  external  fields.  Such  equations  are  of  interest 
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for  two  reasons.  Firstly  for  comparison  with  equations  obtained  by 
other  methods,  and  secondly  because  of  their  close  connection  with  the 
field  equations  for  spinning  media.  The  connection  can  be  seen  in  [3] 
where  we  are  led  to  equation  (1.1)  as  a  consequence  of  Maxwell's  equa¬ 
tions  and  single  particle  equations  of  motion.  The  single  particle 
equations  will  in  fact  be  shown  to  follow  from  the  four  dimensional  action 
integral  of  the  field  theory  by  "variation  of  world  lines". 


It  will  also  be  shown  that  the  Dirac  equation  in  general  rela¬ 
tivity  leads  to  equations  of  motion  and  spin  having  the  same  form  as  the 
Lagrangian  equations. 

It  is  appropriate  at  this  point  to  introduce  some  notation  and 
useful  identities  satisfied  by  a  scalar  Lagrangian.  (For  further  elabora¬ 
tions  see  Appendices  one  and  two.)  We  consider  relative  tensor  fields 

$  (x)  .  Here  A  denotes  both  a  labelling  of  the  tensor  and  its  collection 
A 

of  contravariant  and  covariant  indices.  Thus, 


a,  . .  .a 

{*AW>  s  <*(.)»  £ 


a^. . .a 


$A(x)  denotes  the  same  tensor  with  its  covariant  indices  raised  and  contra¬ 
variant  indices  lowered.  A  repeated  capital  Latin  index,  such  as  in 
$A  $  ,  denotes  summation  over  both  the  tensor  indices  and  the  set  of  ten- 
sors  to  which  A  refers.  Sometimes  an  index  is  printed  as  A  whenever 
A  and  A  denote  different  sets  of  tensors.  For  example,  may  denote 

/v  ri 

both  a  set  of  tensor  fields  and  the  set  of  covariant  derivatives 
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Under  a  co-ordinate  transformation  xP  =  xP(x^)  a  relative  tensor 
transforms  as 


Vx)  =  AaB<>  Vx)  - 


(1.2) 


g 

(A  is  set  equal  to  zero  whenever  A  and  B  refer  to  different  ten- 
A 

sors,  so  that  the  summation  becomes  a  summation  over  the  indices  of  a 
single  tensor  ).  Define  infinitesimal  generators  by 

LX 
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(I  B)  0  E  (XP  =  6P)  . 
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The  (I  B)  °  ’  s  may  be  constructed  explicitly  from  the  recur- 

xx  K 

sion  formulas  (see  Appendix  one) 


(I)  °  =  -w  <5°  for  a  relative  scalar  $  of  weight  w 
P  P 


(i.e. 

<f>(x)  =  <Kx)  1  9x/Bx|W  )  , 

(1.4) 

B$  a 
•  ;P 

(I  B)  a  63  -  SB  6a  , 

K  A  Jp  a  A  a  p 

(1.5) 

aB  a 
A6;p 

Bx  a  ra  ,  <3  r-a  xcr 
(IA  >p  S6  +  6A  6P  *f}  ’ 

(1.6) 

(where  6B  =  1  if  A  and  B  both  denote  the  same  tensor  and  the  same 
tensor  indices,  and  is  zero  otherwise).  Useful  formulas  are 


$  ,  =3  <*>  +  rp  (i.B)  0  <i> 

At  t  A  ax  A  p  B 


(1.7) 


<j>  -  $  =  -RP  (I.B)  0  $_ 

A  yv  A|vy  ayv  A  yp  B 


(1.8) 
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and 


y*>  -  Vx)  -  Op0  WP) 


(1.9) 


for  the  change  of  0.  under  infinitesimal  co-ordinate  transformation 
XP  =  xP  +  £P(x)  .  It  follows  from  (1.9)  that  the  condition  that  the 
function  L(0  )  be  a  relative  scalar  of  weight  w  is  (see  Appendix  two) 


9l 


ir-r—  (I  ^)  °  i  +  W  5"  L  =  0 
8$.  A  p  B  P  ~ 


(1.10) 


Specialising  (1.10)  to  the  case  of  a  scalar  density  (w=l)  L(0A)  ~ 
L($A5^A|a^  depending  on  a  set  of  relative  tensor  fields  $A  and  their 
first  covariant  derivatives  O^i  >  defining  the  variational  derivative 


6l  .  . 

~  _  A  _  Aa 
5$  ~  ~  la 


(1.11) 


where 


A  _ 


9L 


9$, 


and  L 


Aa  _ 


9L 


90 


A  a 


(1.12) 


and  defining  tensor  densities 


a  _  ~a  .  tAQ 

t  =l6  -  0  ■  L 

~p  ~  p  A  p  ~ 


uTa  =  LAT(I  B)  a  o 
«  p  ~  UA  ;P  B 


(1.13) 


gives  the  identity 


6  L 
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(1.14) 


. 


CHAPTER  TWO 


EQUATIONS  OF  MOTION  FOR  A  SINGLE  PARTICLE 


§2.1  Introduction. 

We  will  be  concerned  with  the  derivation  of  equations  of  motion 
for  the  momentum  P  and  spin  S  of  a  single  particle  moving  in  given 
external  fields.  In  special  relativity  theory,  the  natural  definition  of 
an  angular  momentum  tensor  is 


V  *11  Pv  %  py  * 

where  P^  is  the  particle  four  momentum.  is  antisymmetric  and  in  a 

system  of  zero  three  momentum,  i.e.  P^  =  =  ?3  =  ®  »  ^as  on-*-y  three 

non-vanishing  components. 


M 


(0) 

*yv 


o 

0 

0 

-x  P(0) 
T4 


0 

0 

0 

-x  P(0) 
x2  4 


0 

0 

0 

-x  P(0) 
3  4 


x  p(0) 

T4 

x  P(0) 
X2  4 

x  P(0) 

x3  4 

0 


This  may  be  written  as 


=  0  . 


(2.1.1) 
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We  will  consider  particles  which  also  possess  an  intrinsic,  or 
internal,  angular  momentum,  independent  of  a  choice  of  origin,  and  this 
"spin"  will  also  be  represented  by  an  antisymmetrical  tensor,  say. 

Like  the  orbital  angular  momentum,  we  will  assume  that  the  spin  is  char¬ 
acterised  by  three  components  in  a  system  of  zero  three  momentum,  how¬ 
ever  we  this  time  take  the  space  components  to  be  non-vanishing. 


S(0)  = 

yv 


0 

S(°) 

s(0) 

12 

13 

(0) 

0 

s(0) 

12 

23 

(0) 

-S(°} 

0 

13 

23 

0 

0 

0 

0 

0 

0 

0 


i.  e. 

s(0)  pv(0) 
yv 


=  0  when 

y4 


px  =  p2  =  p3  =  o 


r4  *  o  • 


i.  e. 


o  -  s<?>  p4(0) 

U4 


=  S  P 

yv 


So  this  condition  is  written  as 


S  PV  =  0 

yv 


(2.1.2) 


Whenever  the  four  velocity  u  is  parallel  to  the  four  momen¬ 
tum,  (2.1.2)  may  be  replaced  by  the  condition 


S  u  =  0  .  (2.1.3) 

yv 

However,  it  is  usually  assumed  nowadays  that  when  a  particle  has  intrinsic 
spin,  that  the  four  momentum  is  no  longer  parallel  to  the  four  velocity. 
This  was  first  suggested  by  Weyssenhoff  and  Raabe. 

With  P*  no  longer  parallel  to  u^  ,  the  conditions  S  ^  u^  = 

and  S  .  P^  =  0  have  different  meanings.  For  an  extended  body  (see  [4]) 

yX 


0 
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these  are  defining  equations  for  the  world  line  of  the  centre  of  mass. 

However,  whilst  S  ,  =  0  determines  a  single  world  line  which  under 

yA 

X 

certain  assumptions  will  be  a  time-like  geodesic,  u  -  0  does  not 

determine  a  unique  world  line.  We  therefore  take  condition  (2.1.2)  to 

X  _  X 

be  valid,  this  equation  determining  the  velocity  u  in  terms  of  P 


and  S 


yv 


Different  methods  have  been  used  to  derive  equations  of  motion 
for  spinning  particles.  There  is  the  multipole  formalism  for  an  extended 
body  which  consists  of  integration  of  the  conservation  identities 
0  =  T°^  | 0  for  the  total  energy  tensor  over  space-like  sections  of  the 

|p 

particle  world  tube  [5].  Equations  obtained  from  a  differential  geomet¬ 
ric  viewpoint  have  been  found  by  Hans  Kiinzle  [6],  and  there  are  generali¬ 
sations  of  rest-frame  equations  to  arbitrary  frames  [7].  Lagrangian 
equations  have  either  been  obtained  by  a  special  choice  of  Lagrangian  or 
only  in  the  approximation  of  special  relativity  [8].  Before  the  derivation 
of  covariant  Lagrangian  equations  it  will  be  useful  to  discuss  the  general 
form  of  a  Lagrangian  describing  a  spinning  particle. 


A  particle  with  intrinsic  magnetic  moment  y  and  intrinsic 

electric  dipole  moment  d  ,  moving  in  an  external  electromagnetic  field, 

has  an  additional  energy  y*H  —  d*E  where  H  =  magnetic  field,  E  - 

electric  field.  We  have  y H  -  d-B  =  \  MUV  F  where  the  four  tensor 

M^V  is  formed  from  y  and  d  in  the  same  way  that  the  electromagnetic 

field  tensor  F  is  formed  from  H  and  E  .  So,  in  terms  of  a  Lagran- 

y\)  ~  ~ 

gian  L  ,  the  magnetic  moment  tensor  for  a  particle  would  be  defined  as 
MyV  =  2  9L/9F  .  We  will  also  assume  that  the  four  vector  potential 


.. 
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A^  appears  in  the  Lagrangian  only  via  a  coupling  to  the  electric  charge 
current  J  ,  so  that  J  =  3l/3A^  . 

To  give  the  equations  for  both  the  trajectory  and  for  the  spin, 
the  Lagrangian  must  also  depend  on  "co-ordinates"  or  "internal  degrees 
of  freedom"  corresponding  to  spin,  whose  variation  leads  to  the  spin 
equations,  just  as  variation  of  the  space-time  co-ordinates  leads  to  the 
equations  of  motion. 


The  spin  of  a  particle  may  be  described  by  the  rotation  of  a 

set  of  axes  attached  to  the  particle,  so  we  introduce  an  orthonormal 

(a)  (a) 

tetrad  defined  at  each  point  of  the  particle  world  line,  e^  =  e^  (s)  , 

s  =  proper  time  along  the  world  line,  (a)  goes  from  1  to  4, 

e (a)  =  qab  ,  qab  =  (1,1, 1,-1)  .  These  four  four-vectors  constitute 

l" 

the  "spin  co-ordinates"  and  we  consider  a  Lagrangian  that  is  a  function 
of  u^  =  dx^/ds  ,  e^  ,  e^  =  6  e^/^s  (absolute  derivative),  and 
external  fields 


,  A  (a)  ‘(a)  ,  , 

L  -  L (u  ,ey  ,ey  ,\)  • 


The  natural  definition  of  spin  in  terms  of  L  is 

Suv  E  2  e(a)[U  -?TT  •  (2.1.4) 

9ev] 

This  is  the  generalisation  to  relativity  of  spin  in  classical  mechanics 
(see  Appendix  three.) 
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§2.2  Constraints. 


In  the  derivation  of  equations  we  cannot  vary  all  the  co-ordi¬ 
nates  freely  since  they  are  not  all  independent  and  must  satisfy  con¬ 


straints  on  variation.  To  obtain  the  spin  equations  we  vary  the  vectors 

(a) 

e^  keeping  them  orthonormal,  i.e.,  we  vary  e^  subject  to 

e^a^  .  The  variations  de^  are  not  all  independent  and 


_  ((a)  (b))y 

must  satisfy  de .  e 


(a) 


=  0  . 


'y  -  =  0  ,  or  equivalently  de^  e(a)v) 

The  constraints  are  ten  relations  between  sixteen  variables  so  we  have 


essentially  six  independent  variables. 


Introducing  scalar  Lagrange  multipliers  y^ (s)  =  y^a(s) 


and  L  =  L  +  y^Ce^3)  e^^  -  ,  then  extremising 


L  ds  for 


arbitrary  variations  de 


(a) 

y 


vanishing  at  the  endpoints  gives 


SL 


6e 


(a) 


=  0 


i.e. 


6e 


=  0 


y  y 

The  Lagrange  multipliers  are  therefore 


(2.2.1) 


y 


1  6L 


ac  2  (a)  (c)y 

oe 

y 


To  eliminate  y  ,  in  (2.2.1)  and  obtain  spin  equations,  multiplying 

ab 


(2.2.1)  by  e 


(a) 


v 


and  antisymmetrising  gives 


6L  (a)  _  9  ,_(b)y  (a)  (b)  (a)y. 

e  ,  =  -2  y  .  (e  e  .  -  e  .  e  ) 


6e 


(a)  v] 

[y 


ab 


V 


V 


(b)y  (a)  9  (a)y  (b) 

=  -2  u  ,  e  e  +  2  y  ,  e  e 

Mab  v  ab  v 


=  0  (since  uab  =  yfca  )  • 


’ 
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These  are  the  six  spin  dynamical  equations. 


Variation  of  the  world  line  xy  (s)  ,  s  =  proper  time,  to  obtain 


equations  of  motion  does  not  lead  simply  to  6l/6x^  =  0  where  the  action 


is  I  = 


L(xy  ,  dx^/ds)  ds  because  of  the  constraint 


^  '  -1 


* 

(x*"1  =  dx^/ds)  .  We  remove  this  constraint  by  a  process  of  "parametrisa- 
tion".  Introducing  an  arbitrary  invariant  parameter  A  ,  s  =  s(A)  ,  then 


f  ,  y ,  v  1/2  (  1/2  , 

ds  =  (-SyVdX  dX  )  =  X  ^  ^ 


f  •  dX 


•u 

where  x 


=  dx^ /dA  .  We  rewrite  I  as 


I  = 


L(x^*  (s)  jdx^/ds)  ds 


L(x^(A),x^  dA/ds) 


dA 


L(x^,x^,ds/dA)  dA 


defining  a  new  Lagrangian  L(x  ,x  ,ds/dA)  .  The  variables  of  L  are 

*  u  • V  ds  2 

still  constrained  since  we  now  have  g  x  x  =  replacing 


d/  dxv 


g  i“L-  =  -1  .  We  eliminate  the  constraint  by  using  the  above  equa- 

Byv  ds  ds 

r*j 

tion  to  replace  ds/dA  in  L  .  Then 


I  = 


L(xy,xy,(-gyvxMxv)J’/")  dA 


•yv,l/2, 


L1(xy,xy)  dA 


where  is  defined  as 


. 
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y  *y  ~  y  *y  *yv.l/2. 

Lx(x  ,x  )  =  L (x  , x  , ("gyvx  x  )  ) 


,  *yv.  1/2  ,  y  ,  #y*Vv-l/2  *y. 

=  (~8yvX  x  )  L(x  »("gyvX  X  ^  X  ) 


•y 

The  Lagrangian  is  homogeneous  of  degree  one  in  x  ,  irrespective 

of  the  dependence  of  L(x^,dx^/ds)  upon  dx^/ds  . 


Our  action  principle  will  henceforth  be 


0  =  6l  =  6 


L^(x^,x^)dA 


for  arbitrary  variations  in  x^  (no  constraint)  where  A  is  any  invar¬ 


iant  parameter  and  is  homogeneous  of  degree  one  in  x 


;y 


§2.3  Covariant  Equations  of  Motion  in  Given  External  Fields. 


61  =  6 


The  equations  are  obtained  from  the  action  principle 

dA  =  0  for  arbitrary  variations  of  x^(A)  and  e^(A) 


with  fixed  endpoints  x^(A^)  »  e^Q^)  ,  (i  =  1,2)  . 

r* 

iant  parameter,  \  =  X(s)  ,  and  is  homogeneous  of  degree  one 

v^  =  dx^/dA  . 


\  is  any  invar- 


in 


Ln  =  L  (vy,e^a\6e^a^/6s,$A) 

1  1  y  y  o 


(2.3.1) 


The  set  consists  of  the  external  fields  <p^  ,  the  Riemann  tensor 


.a 


3y6 


and  their  covariant  derivatives. 


$A  ’  R01gy6  ’  |  a  ’  r016y6  |  u 


,  •  •  •  )  • 


. 


a£o 


■ 


' 
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Defining  e^  =  /6\  ,  then 


6e(a)/6s  =  4(a)  dX/ds 

y  y 


•(a)  ,  (b)  (c)  a  3. -1/2 

e  (-TL  e  e^  v  v  ) 

y  be  a  3 


so  may  be  reexpressed  as  a  new  function 


L,  =  L(vW,e(a),4(a),$A) 

1  y  *  y  A 


We  define  the  canonical  momentum  P  ,  the  spin  angular  momentum 

y 

Spa  =  SCpa]  and  the  multipole  moments  by 


P  =  ,  s  =  2  efa>  -4rr  .  ^  -  fr 

p  [p  9l(a)  9<*\ 


y  3v^ 


'a] 


(the  derivatives  3L/3$a  =  M~  are  defined  such  that  dL  =  M~ 


d$> 


has  the  same  algebraic  symmetries  as  .  For  example,  if 

L  =  F 


XPV  with  F  antisymmetric,  then  L  =  F  .  X^PV^  and 

yv  yv  J  yv 


is  defined  as  ). 


In  terms  of  L-^  we  have 


3$a 


(a) 


,  S  a  =  2  er 

P  [p 


8L1  .  dA 

(a)  ds 

3  — 2 — L 


3L^  3L^ 

Py  "  .  u  +  r^(a)  v  6s 


6s 

(a) 


9  (^— ) 


3v 


6e  3vh 

q 

6s 


8L1  .  8L  1  • (a)  ,dXx  3 

+  - 7-^r  e„  v.,  (— ) 


3vP  6e 


(a)  ca  vy  Ms' 


3 


a 


6s 


Since  is  homogeneous  degree  one  in  vP  we  have  S  °  and  P 


y 


(2.3.2) 


(2.3.3) 

and 

3L/3F  . 

q3 


(2.3.4) 


(2.3.5) 


both 


. 


.  , 


i  ' 
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independent  of  parametrisation,  and  proportional  to  ds/dA  . 

(  3. ) 

A  variation  de /j  in  the  tetrad,  keeping  the  world  line 

p 

fixed,  gives  (since  g^v  and  are  unaltered)  the  variation  of 

*  (a) 


ey  '  to  be 


*  (a)  ,  (a) 

w  y  x  _  wdey  r  3  (a)  An  _  6  ^  (a), 

d(  SA  }  d(  dA  ry  A  eB  v  )  5A  (dey  }  * 


Hence 


dL  =  de.(a)  +  3L 


(a) 

;y 


(de,,  ) 


(a)  6A  wcy 


'y 


,_9L  j5_  JL  ..  (a)  ,  _6_  9L  (a). 

”  (.  (a)  “  6A  ^  •  (a)  }  dey  +  6A  (,-(a)  dey  } 
de  de..  de. 


U 


11 


11 


9L  ,  (a)  . 

- 7-7-  de  is  a  scalar,  so 

y 

y 

spin  equations  are  therefore 


6  8L  ,  (a).  d  (  dL  i  (a)\ 

(~ T77V  dei,  )  =  7T  (“TTIT  de,,  )  . 


6A  (a)  y 

3eu 


dA  v(a)  y 

3ey 


The 


5l  (a) 


(a)  '  V] 

[y 


=  0  where 


6l  9l 


(a)  .  (a) 


_5_  _8L 

<5A  (  • 


'y 


'y 


9e. 


(a) 

y 


)  .  (2.3.6) 


Identity  (1.10)  with  w  =  0  for  the  particle  Lagrangian  is 


P/  -  +  /-(I.®)  V  *.  -  0  •  (2.3.7) 


A  '\i  ‘B 


Rearranging  and  antisymmetrising,  this  is 


p'v1 


-  e 


(a)  [y  <$L 
6e  V] 


_6_  ,  (a)[y 

6A  v 


9e 


9L  j.  mA/t  B^  tyv], 
TS)5  +  m^(ia~)  *b 


=  0 


v] 


Spin  equations  (2.3.6),  together  with  definition  (2.3.3)  for  and 


the  above  identity  then  give 


15. 


(2.3.8) 


To  obtain  momentum  equations  we  make  an  arbitrary  variation  in 
the  world  line  dx^(A)  (dx^  =0  at  endpoints)  keeping  e^a^  covariantly 


y 


(a) 


constant  on  variation,  6ev  '  (A)  =  0  .  The  most  elegant  approach  is  a 


y 


manifestly  covariant  one. 


The  change  in  L  is  dL  =  6L 


3L  yN  ,  3L  ^_(a)  .  3L  j 


3v 


y 


6(vM)  + 


3e 


(a)  y 


6eva/  +  ■WJ/  \  6(evcl/)  +  X't*  <5$ 


'y 


3e. 


(a)  y 


'y 


3$  A 
A 


where  6 (vy)  etc.  are  the  absolute  changes  in  each  dynamical  variable, 


,dx 


y 


,dx 


y 


v 


6(VH)  -  6<^)  -  d(^)  +  r/Y  %  dx^  =  ±  (dxU)  . 


*  (a)  -  n 
oe  =0  , 

y 


5  $>  =  $A ,  dxy  . 

A  Ay 


Owl  Owl 


,  ,  6e(a) 


-R  e<a)Y%  dx° 
yypa  dA 


3L  <w  ys  6  ,  3L  ,  Us  6  ,  3Ls  j  u 

6P 

-jr-  (P  dxy)  -  — r^-  dxy  since  P  dxy  is  a  scalar 

dA  U  6A  y 


dL  dA  = 


6P 


(- 


J±  


6A 


R  ft  - --Vy  vY  +  M~  $Ai )dxP  dA 

a3yy  -sl(a)  A  |  y' 


3e 


a 


a 


By  definition  (2.3.3)  for  S  ,  the  momentum  equations  are  therefore 


»3f  {'3  /lJ-pqhfIS  3fi 


. 
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— kL  =  I  r  s°^  $  | 

6A  2  Vyy  A  y 


(2.3.9) 


In  the  first  term  we  see  the  interaction  between  the  spin  and  the  gravi¬ 
tational  field. 


Together  with  (2.3.8)  and  (2.3.9)  we  need  an  equation  determin¬ 
ing  which  we  take  to  be  =  0  .  In  the  simplest  case  with 

=  {0}  ,  using  (2.3.8)  and  (2.3.9)  gives 

0  =  £  (S“P  V  ■  2p[V5Ipe  +  s“e  i  Vpm  vY  sW) 


i.e.  Pa(v%)  -  va(P%)  -  \  SaB  R6ypu  v 


Y  SPU  =  0 


(2.3.10) 


s  „  s 

yv 


yv 


is  conserved  since 


£  (S  ,SUV)  -  2  S  % 

6A  v  yv  yv  6A 


yv 


,  y  v  v  U, 
2  S  (P  v  -  P  vM) 
yv 


=  0 


§2.4  Motion  in  an  External  Electromagnetic  Field, 


As  a  special  case  of  §2.3,  take  -  {A^,A  where  R 


B 


,a 


denotes  R~  r  and  A  is  the  four  vector  potential.  Assume  that  A  . 

8yo  a  a  I  p 


appears  in  the  Lagrangian  only  through  F^  =  A^ | ^  -  A^ | ^ 
the  magnetic  moment  as 


We  define 


m“B  = 


9L 


9A 


=  2 


9L 


$  a 


9F 


a$ 


■  1 1 1 5  ■■  ■  •  - 


■ 
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A  $Y<5 

and  the  quadrupole  gravitational  moment  tensor  as  Q  =  Q 

LI 

(X 

3L/3R  £  ^  and  we  assume  the  dependence  of  L  upon  is  given  by 

3L/3Aa  =  ja  =  e  va  .  (e  =  electric  charge) . 


M^(I  °  $ 

;P  B 

~  (-<$a66)A R  +  -  (-6^16a  6^2  -  6a  63l6^2)Aft  ,r  +  QA(I  B)  V 

3Aa  *\\*2  ot1  p  p  ^  A  p  *B 


3l  a  b  a 

BAa  ~p  "  3Aa  la  ~allp  "  8Aa|a.  pla2  A  p  B 


8L  A  . 


a  aa  3y<5  a  „  ^  oy&  a 

-evA-M  F  +  Q  1  Ro.-30  R  x 

p  pa  xp  3y<5  a  PY<5 


...  =  -e  A[Pv01  -  F[P  M°]a  -  4  RtP_,  Qala6Y 


B 


a 


a3Y 


Substitution  of  the  above  into  (2.3.8)  and  defining  a  new  "kinetic"  mom¬ 
entum  p  =  P  -  e  A  ,  gives  spin  equations 
y  P  P 


(2.4.1) 


6P 

6X 


- 


$ 


Ay 


6Py  V 

— xy  +  e  A  |  v  — 

oA  y  v 


. a  i  -  — — —  Aii  -  QAR  i 

3AV  Av|y  3Aa .  aA | v | y  p  A|y 


6P 


^  +  e  F  vv  -  ~  Mv*  F 


6A  vy 


_  q  3y5  Ra 

vX  y  a  3y<$  I  y 


The  momentum  equations  (2.3.9)  are  therefore 


“5x  =  I  Ra6yw  SaBvY  +  e  VV  +  7  mVXfvX|u  +  Qa6Y6  R°M|y 


(2.4.2) 


'  i-  '  .  -  : 


i  3  t  - 
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In  the  absence  of  anx  fields,  =  0  =  ^  =  dX  ’ 

(2.4.1)  and  (2.4.2)  reduce  to  dS^V/dX  =  p^vV  -  p^v^*  and  dp^/dX  =  0  , 
and  (2.3.10)  gives  p^  parallel  to  v^  .  This  implies  dS  /dX  =  0 

and  dJ^V/dX  =  0  where  J^V  (orbital  angular  momentum)  is  = 

y  V  V  y 

x  p  -  x  p 


.  .  y  , 

Now  choosing  X  to  be  proper  time,  X  =  s  ,  v  v  -  -I  ,  we 
can  multiply  (2.4.1)  by  and  solve  for  p^  (set  =  0  for 


a 


simplicity) .  We  obtain 


y  y  .  [y  v]a  6S^V 

p^  =  mvM-2FL  M  v  -  —7 —  v 

F  o  a  v  os  V 


(2.4.3) 


(mQ  =  -p^  v^)  .  Notice  that  p^  is  not  parallel  to  v^  ,  p^  is 

y 

expressed  in  the  above  equation  as  a  sum  of  a  term  parallel  to  v  (i.e. 
m  v^  )  and  a  term  orthogonal  to  v1"*  .  A  satisfactory  feature  is  that  the 


o 


energy  momentum  vector  p^  should  also  depend  on  the  spin. 


Now  replacing  p^  by  the  right  hand  side  of  (2.4.3)  in  (2.4.1)  gives 


,yv 


=  -2  v  -------  vV]  -  2v  F^  M^°V^  +  2v  F^  M^avV-*  -  2F^y  MV^a 

6s  X6s  Xa  Xa  a 


i.e. 


_  2v  v]  .  2v  MXaF  [yvv]  +  2v  pAc^tu  „vl  _  2FIy  Mvla 

X  a  a 


:yv 

r  -  2\ 

6s  X  6  s 


X  a 


In  the  case  M°^  =  K(s)  S°^  ,  =  0  ,  we  have  the  spin  con¬ 


servation  equation  -^7  (S^S01^)  =  0  from  (2.4.1). 


' 


■ 
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In  summary,  sections  (2.3)  and  (2.4)  contain  covariant  Lagran- 
gian  equations  of  motion  for  spinning  particles  in  gravitational  and  any 
other  external  fields.  The  ease  of  derivation  and  their  generality  are 
features  of  equations  (2.3.8)  and  (2.3.9).  They  are  general  in  the  sense 
that  no  explicit  dependence  upon  its  variables  is  assumed  for  the  Lagran- 
gian,  L  is  any  function  of  v^  ,  e^a^  etc.  The  only  possibility  for 

r™ 

obtaining  different  equations  might  be  by  use  of  different  definitions  of 
spin,  four  momentum  etc.  (though  the  content  of  the  equations  would  not 
change).  For  example  one  may  wish  to  define  a  new  momentum  in  (2.4.1) 
and  (2.4.2)  by  P^  =  p^  +  q^  where  q^  is  the  electric  dipole 

moment.  (2.4.1)  and  (2.4.2)  then  become  (with  Q  =  0) 

ot 


6p 


y 


6X 


-x-  R  S0^  vY  +  e  F  vV  +  \  MVA  F  .  .  +  (F  .  qX) 

2  yya$  yv  2  vA|y  6A  v  yA  H 


^  =  2P[V]  -  2F[lJ  (MV]“  +  vVl  q“)  . 

0  A  Ot 


These  are  the  equations  obtained  by  Suttorp  and  de  Groot  [9]  in 
the  special  relativistic  case. 


. 


■ 


■ 


CHAPTER  THREE 


RELATIVISTIC  FIELD  EQUATIONS  FOR  POLARIZED  MEDIA 


§3.1  We  consider  a  polarised  dustlike  medium  interacting  with  the  grav¬ 
itational  field  and  any  other  fields  present.  The  medium  is  taken  to 
consist  of  coherently  moving,  spinning  particles  which  interact  only 

through  the  external  fields.  The  particle  world  lines  are  the  integral 

cx  a 

curves  of  a  numerical  flux  vector  field  n  u  (x)  ,  u  u  =  -1  ,  describing 

UG 

a 

the  medium,  which  is  assumed  to  satisfy  (nu  )i  =  0  (conservation  of 

I  Uu 

number) .  The  internal  spin  of  the  medium  is  described  by  an  orthonormal 

(a) 

tetrad  field  ev  '  (x)  which  also  serves  to  specify  the  metric  field. 

a 

Variation  of  a  four  dimensional  action  integral  leads  to  gravi¬ 
tational  equations  and  equations  for  the  spin,  by  variation  of  the  tetrad 

fa') 

field  eK  '  (x)  .  The  resultant  spin  (field)  equations  are  equivalent  to 

CL 

the  single  particle  spin  equations  of  motion,  and  by  a  variation  of  world 
lines,  the  four  dimensional  action  principle  also  gives  the  single  particle 
translational  equations  of  motion.  Non-gravitational  fields  present  will 
of  course  have  field  equations  obtained  by  variation  w.r.  to  their  compo¬ 
nents.  The  outcome  is  thus  a  complete  set  of  field  equations  relating  the 
fields  to  their  sources  and  also  the  single  particle  equations  obtained  in 
chapter  two. 

The  equations  are  obtained  from  the  action  principle  6l  =  0 
for  arbitrary  variations  in  the  dynamical  variables,  I  being  a  four 
dimensional  action  integral. 


20  - 


■ 

. 


■ 


■ 


■ 


. 


I 


(3.1.1) 


(-C16TT)  1(-g)1/2R  +  L($A,$A|a))d4x 

rsj 

-1  1/2 

We  have  taken  —  (16tt )  (-g)  R  ,  where  R  is  the  curvature  scalar,  for 
the  gravitational  free  field  Lagrangian.  The  scalar  density  L  is  a 
function  of  the  fields 


$A  "  {2a’eaa)’VRa8Y«} 


(3.1.2) 


Ct  ( 3. ) 

where  N  =  /-g  nu  is  the  numerical  flux,  e\,  is  the  tetrad  field, 

^  °  CL 


a 

R  ^  ^  is  the  Riemann  tensor  and  ip^  an  arbitrary  set  of  fields.  We  are 

assuming,  for  simplicity,  that  L  does  not  depend  on  second  and  higher 

covariant  derivatives  of  the  $A  (the  general  case  with  dependence  on 

higher  derivatives  is  discussed  in  Chapter  Four).  We  also  assume  that  no 

Oi  OL 

derivatives  of  R  ^  ^  and  N  appear  in  L  . 


Before  proceeding,  a  few  comments  on  the  role  of  the  tetrad 

C  3  ) 

field  is  in  order.  Firstly,  the  components  e^  (x)  are  the  gravitational 
variables  just  as  any  orthonormal  tetrad  field  describes  the  metric 
g  =  D  ,  e^e^  .  Secondly,  as  a  set  of  spin  axes  for  each  world  line 
they  are  spin  variables  that  will  be  dynamically  determined  along  each 
world  line  in  accordance  with  the  particle  spin  equations  of  motion.  Be¬ 
cause  of  this  second  role  the  Lagrangian  has  an  explicit  dependence  upon 
the  in  addition  to  its  dependence  via  the  metric  g^  (so  that  we 

could  not  replace  e^3^  by  8^  in  (3.1.1)  and  (3.1.2))  and  hence  L 
would  not  be  expected  to  be  invariant  under  an  arbitrary  field  of  tetrad 
rotations.  One  could  have  considered  a  Lagrangian  density  depending  upon 
both  g  and  e^a^  (and  then  varied  g  to  obtain  the  ten  gravitation- 
al  equations  and  varied  six  independent  components  of  the  tetrad  to  obtain 


■ 


•  .x.,!  }  ftiaiavoo 


“ 
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spin  equations).  However,  since  the  orthonormality  conditions  g  = 

r  V 

r|  ,  e^a^e^^  are  ten  relations  between  the  twenty  six  variables  g,„,  , 
ab  y  v  yv 

(a) 

e^  ,  we  have  sixteen  independent  variables  and  it  is  much  simpler  to 
(cl) 

take  the  e^  as  our  sixteen.  We  take,  therefore,  arbitrary  independ- 

ent  variations  in  the  e  to  obtain  ten  gravitational  equations  and 

OL 


six  spin  equations. 


Much  guidance  for  the  following  calculation  is  due  to  Rosenfeld 
(see  [10])  (and  also  Belinfante  [11]).  The  difference  between  their  work 
and  the  present  is  that  spin  is  considered  in  their  scheme  only  in  the 
sense  that  the  fields  are  allowed  to  be  spinor  fields. 


§3.2  Tetrad  Variation. 


We  make  an  arbitrary  variation  de 

(a) 


resultant  variations  in  functions  of  e 


y 


(a) 

y 

are 


details) 


in  the  tetrad  field.  The 
(see  Appendix  four  for 


dg  n  =  2  T|  ,  e[a^  de^!  for  the  metric  , 
&a$  ab  (a  3) 


(3.2.1) 


d  r*3  =  dgP  N  —  -r-  dg__  for  the  Christoffel  symbols,  (3.2.2) 

ax  (a|x)  1  CTT 


dRa^,  c-  =  2  dr^rr.i  ,  for  the  Riemann  tensor 

3Y<$  3[6|y] 


(3.2.3) 


It  follows  from  (3.2.3)  that  for  an  arbitrary  tensor  density  A 


3y<$ 


a 


that 


■ 


noisi.  u.\u'i.  I  'c  woi lo\  aria  30"-  &oi:K>biap  ti:  if 


- 

V 


23. 


By6  a 

A  PY  dR  0  x  =  -2  A 
a  8yo  a 


B[y6]  ra  fA  . 
i  dr  OJC  +  (div) 
y  Bo 


(3.2.4) 


where  (div)  represents  a  divergence  term  (which  by  Gauss’s  theorem, 


will  play  no  part  in  the  derivation).  From  (3.2.2)  it  follows  that  for 


TO 


an  arbitrary  tensor  density  B  that 


drPCTT  =  \  (Bp^  +  Ba[pT]  -  BT(P0))|t  dgpa  +  (div)  .  (3.2.5) 


The  variation  in  the  free  field  gravitational  Lagrangian  is 


(div)  +  (— 16TT )  1  d((-g)1/2R)  =  (lev-1  (-g) 1/2  Gae  dgaB 


=  (87T)  1(-g)1/2  GaS  e(a)ade^a)  .  (3.2.6) 


Care  must  be  taken  in  calculating  the  change  dL  in  L  since  L  depends 
,  (a)  . 

upon  the  e^  in  three  ways : 


(i)  L  has  an  explicit  dependence  on  e 


(a)  ,  (a) 

u  and  eu|v  ; 


(ii)  L  depends  on  R 


a 


By  <5  ’ 


(iii)  L  depends  on  the  T  ^  's  hidden  in  covariant  derivatives 
0.  I  (i.e.  e^a]Q  and  ip A  i  ). 


A  a 


a  B 


A  a 


,  %  ^ 

The  change  in  L(4>,$.i  )  due  to  (iii)  is  (3L/90  i  )  ~ 

~  A  A  CX  Ot  r\  -p  P 

~  1  a  1 


rv  1 


Om  JT,p 


OT 


which  by  (1.7)  is  (9L/30a  ,  )  (I  ~)  drM  .  Defining 
]  ~  A  t  A  p  B  OT 

I  r*j  1  nu 


dr 


OT 


L~  =  9L/90.  ,  L~a=aL/3$., 

~  ~  A  ~  ~  A  a 


6l 


60 


A 


=  lA  -  lA*  ,  UT°  =  lAT(I  *)  °  *_ 
a  ~  p  A  p  B 


■ 
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3y6  a 

as  in  Chapter  One,  and  defining  =  9^/3 R  ,  then  the  change  in 


3y<5 


L  is 


6l 


(a)  ,  tttg  *Pp  ,  3y6  Jna 


9L 


dL  =  — dew  +  Ulun  fir*'  +  0 
~  -(a)  y  ~  n  gt 


p  gt 


dR  By6  +  (,  (a)  “'a  ' |a 

de 


de^a)V  (3.2.7) 


y 

6l 


g  a 


-  *  s,sl,i'i,  <■“«  *  «*•> 


(3.2.8) 


O  £ 

using  (3.2.4)  and  (3.2.5)  .  Antisymmetry  of  £a  Y  in  the  first  two  and 
last  two  indices  and  another  use  of  (3.2.5)  gives 


qBy6.  dr“  =  -q  T(po)Y,  dg+(div)  . 
2  a  y  36  -5  |yt  6pg 


Using  (3.2.1),  (3.2.6)  and  (3.2.8)  the  action  principle  6l  =  0  for  any 

C  ) 

de  therefore  gives 

y 


—  ( 8tt ) — 1  (-g)1/2  GW  e  -  i 

6ev 


+  (U  P[VT]  +  uv[pT]  -  uT(pv) 


»N/  ^ 


)  i  e  /  s 
'  T  (a)p 


+  42T(PV)Y,  e 


YT  (a) p 


(3.2.9) 


or  equivalently. 


(8tt)  1(-g)1/VV 


6L 


6e 


_  e(a)y  +  (UM[VT]  +uv[pt]  _  DX(WV) 

(a)  ~  ~  ~  |T 


v 


+  4  Q 


x(yv)Y 


YT 


(3.2.10) 


The  above  are  the  gravitational  and  the  spin  equations.  Since  every  term 


. 


' 


. 
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is  symmetric  in  y,V  ,  except  for  e 
of  motion  are 


(a)U 


6l 


(a) 


,  the  six  spin  equations 


V 


(a)  [y 


6L 


(a) 
:  v] 


=  0 


(3.2.11) 


The  next  step  is  to  replace  e 


(a)y 


6L 


6e 


(a) 


in  (3.2.10)  by  terms  having 


V 


possible  physical  interpretation  and  this  is  achieved  by  means  of  the 
identity  of  the  type  (1.14)  that  L  ,  as  a  scalar  density,  must  satisfy. 
Recalling  (1.14)  which  reads 


6l 


TG  G 

U  |  +  t  +  ..  , 

~  p  T  ~p  0$A  A  p  B 


(I  a  <&  =  0 


(3.2.12) 


f  G  _  rG  *  _  A O. 
(t  =  L  6  -  $.  |  L~  ) 
~p  ~  p  A  p  ~ 


(3.2.13) 


in  the  present  context  with  L  =  )  , 


/■v;  rv 


r^j  r+j 


°A  ,ea  ),l^A,R  3y<$^  *  °A|a  ^ea|3*^A|c^ 


(3.2.14) 


we  have 


6L  r  - 

~  /  ~r  Bn  G  ^ 

MTVp  *2 


3L  ^  3L  „  ^ 

Na  -  -2—  Na  6a  - 


6L 


8# 


3N 

6L 


.a 


~  (a) 

—  e 


p  6e(a)  P 
6G 

«aV  \  +  sXV  rb 


6^a  A  p  B 


(3.2.15) 


(R  is  short  for  RaD  ~  ,  and  written  out  explicitly  the  last  term  is 


B 


$Y<5 


.A,,  Bv  g 


2  <XA  >D  *B  “  S. 


Ro 


p  b  ^p  3y<5  py<$ 


-  32, 


gy  6 


sOt 


(3.2.16) 


s>-5»  flOi3«*r*  So 
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The  non-gravitational  field  equations  are 


6L 


Sip, 


=  0 


(3.2.17) 


(3.2.15)  now  gives  (modulo  the  non-gravitational  field  equations  (3.2.17)) 
the  identity  (3.2.12)  in  the  form 


6l 


~  e(a) 


3L  _  3^ 


6e(a)  P 

o 


Replacing 


9N 


P 


N0  - 


8N 


Na6a  +  UTa  |  +  t  Q  +  QA(I  B)  0  R  .  (3.2.18) 

a  ~  p  ~  p  T  ~p  ~  v  A  p  B 


6L 


6e<a) 


e^a^  in  (3.2.10)  by  the  above  expression  gives  the 


gravitational  equations  in  the  form 


(8,rVg)1/2  gpv  =  spu  -  \sa  8UV  +  iyv 

3NP  3N 


+  (UW[VTl  +  UV[WT]  +  uT[vy]) 


+  4gT(yV)Y,T  +  SA(IAB)yV  Rb 


(3.2.19) 


§3.3  Lagrangian  Decomposition. 

We  now  write  L  of  (3.2.14)  as  a  sum  L  =  L..  +  L0  where  L0 

IN'  '■N'  ^  ^  Z.  ~ 

is  the  free  field  Lagrangian  for  ip^  and  represents  the  matter  and 

its  interaction  with  the  applied  fields. 


L  =  L1 (Ny 


>V  +  ^A’4' 


A  a 


) 


(3.3.1) 


' 


Mn  s\r. 


. 


27. 


•(a)  _  ^ (a)  ji 
a 


e  ■  I ,  uyu  -  -1  ,  ip.  =  {ipA  |  ,R^0  J  .  (3.3.2) 

a|y  y  A  rA,rA  a  3y<$ 


.a 


is  assumed  to  be  homogeneous  of  degree  one  in  its  first  argument 

„y  ) —  y 

N  =  n/-g  u  so  that 


r  / —  t  /  h  (a)  •  (a)  .  N 

J*1  =  n/-g  Ll(u  ’ea  *ea  ’V 


(3.3.3) 


where  L^(uy,...)  is  obtained  from  L^(NM,...)  by  simply  replacing 


y 


n1 


y 


y 


with  u  .  This  assumption  will  be  clearly  satisfied  since  for  dustlike 


material  the  flux  of  material  four  momentum,  and  of  spin,  (defined  as 
derivatives  of  L^(Ny,...)  )  will  be  proportional  to  the  particle  density 


n 


is  a  single  particle  Lagrangian  as  in  (2.3.1)  with  parameter  A  set 
equal  to  proper  time  s  .  We  have 


2cM  5  typers  ■  M 


a 


(3.3.4) 


where  Q  is  the  quadrupole  moment  per  particle.  Define  a  "spin  flux" 


syxv  =  2  ^[yx] 


(3.3.5) 


Then,  with  $A|  =  4A|„,e(^}  ,  we  have 

A|a  A|a  a  3 


SMTV  .  „V[|IT]  .  2(3L/34  )(I  5)ttp]  * 

~  ~  ~  A  v  A  B 


=  2(3L/3*A|v)(IAB)[TlJ^B  +  20L/3e^)(-6^g6T])e<a) 


=  s'™  -  20Ll/3ef(a]  )  e(a)Tl 
~(\p)  ~1  [y  I  v 


■ 


V 


' 


. 


.7  ;  - 
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(defining 


yxv  =  v[yx] 
-Of,)  ~  -Of,) 


E  20L/3*A|v)(IAB)[Ty]  *B 


the  spin  flux  of 


the  field) . 


The  spin  flux  therefore  consists  of  a  field  and  a  matter  part, 


sPtv  _  syyv  +  yxv 

—  -(if;)  (mat) 


(3.3.6) 


where 


syTV 

~(mat) 


-2<3V3e[u|v)e 


(a)T] 


o/ /  (a)  \  (a)T]  o/ —  rzr  /c/(a)>t  v  (a)T] 

=  -2/-g  n(dL1/3e^|v)e  =  -2/-g  n(dL1/3e^  )u  e 

-  /~  cT^  V  /—  c^T  V  /o  q  -7\ 

=  -v-g  n  S  u  =  n/-g  S  u  (3.3./; 


UT 

S  being  the  spin  angular  momentum  per  particle  (see  definition  (2.3.4) 
with  X  =  s  ).  Equations  (3.3.4),  (3.3.5)  and  (3.3.6)  describe  the  last 
three  terms  in  the  gravitational  equations  (3.2.19).  All  that  remains  is 
to  examine  the  content  of  the  other  terms  forming  the  right  hand  side  of 


(3.2.19). 


First,  noting  that  L  is  a  function  of  e 


L,  is  expressed  as  a  function  of  e 

-1  a 


(a) 


(a) 

a  I  8 


(see  (3.2.14))  whilst 


since 


•(a)  _  (a) 

e  =  e  i. 

a  a  8 


u 


3 


=  e(^  N3  /  (-N  NA)1/2  , 

a  8  ~  ~*~ 


we  have 


3L  3L- 


3L. 


3r 


3Nh 


+ 


3e 


e(a) 
(a)  a  8 


(-N  NA)  1^2(63  +  u3u 
~A~  p  ( 


a 


(3.3.8) 
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3L  3L 
1  +  L 


(  \  ^  +  U^U  ) 
3up  3e(a)  P  P 


(3. 


a 


TXx  1/2 


(by  the  homogeniety  of  =  (-N^N  )  ) 


^  ot  a 

—  s  = -£ s  =  h 

3N  3n  1 

/V  <■**>» 

3L. 


(3. 


3L  ^  gPH  .  h,  _NV  gPP  +  e.(a)  |w  NV  +  4X  e.(a)uy  NV 


3n 


3uh 


3;<a>  ° 

a 


3e(a)  01 


'a 


3l. 


de 

a 


(3, 


where  PP  is  the  canonical  momentum  per  particle  defined  by  (2.3.5) 


(with  A  set  equal  to  s  ).  tP  defined  by  (3.2.13)  is 


UV 


E  L  gyV  -  $> 'y  3L/3$ 

r+u  A  ~ 


A  V 


vt  ,>r  -  V”  n/»A|V  -  -ru  v».“ 


L1  g  +  <-g) 


e<a>  ly  — At-  NV 
a  (a)  ~ 

de 

a 


(3, 


where  the  canonical  energy  tensor  for  the  fields  ijj.  is  defined  by 


<-8>1/2  'W  “  h  ^  *  V"  ^»*A|v 


.yv  .  y 


(3 


Equations  (3.3.10,  11,  12)  together  give 


dL  3L 

—  NV  gPP  -  — 


3Nh 


3N 


a 


f  gUV  +  tVV  =  PW  NV  +  (-g)1/2 


3.9) 


3.10) 


3.11) 


3.12) 


.3.13) 


.  (3.3.14) 


. 
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Substitution  of  (3.3.4),  (3.3.5)  and  (3.3.14)  into  (3.2.19) 
gives  the  final  form  of  the  gravitational  field  equations  for  polarised 
dustlike  material: 


“1  „yv 


nW  _ 


,y  v 


yv 

'010 


—  (8tt)  -  Gpv  =  Tpv  =  n  PHUV  +  tP:N  +  i(-g)  ^  ‘  (_S/“  V+S"  lp+S'p‘) 


-1/2 ,  lixv ,  0vxy ,  _vyi. 


T 


(3.3.15) 


yv  uv 

T  ,  coupled  above  to  the  Einstein  tensor  G  ,  is  the  symmetric,  covar- 

iantly  constant,  total  energy  momentum  tensor  for  material  and  fields.  It 

is  written  as  a  sum  of  a  material  energy  momentum  tensor  (i.e.  a  momentum 

flux)  n  P^  uV  ,  a  canonical  energy  tensor  t^^  for  the  fields 

defined  by  (3.3.13),  a  spin  term  where  the  spin  flux  S^T^  is  a  sum  of  a 

field  and  a  matter  part  S^TV  =  +  (-g)1^  n  S^T  uV  (see  3.3.6,  7) 

and  a  contribution  from  gravitational  quadrupole  moments. 


§3.4  Single  Particle  Equations  From  the  Four  Dimensional  Action  Integral. 


(a)  [y 


6e 


It  is  easy  to  show  that  the  spin  equations  of  motion  (3.2.11),  i.e, 

6l 

=  0  ,  are  equivalent  to  the  single  particle  spin  equations 


(a) 

v] 


(a) 


(2.3.6).  First,  since  depends  on  e^  only  via  the  metric. 


6L, 


3L, 


3L, 


2e 


3L, 


2e 


6e 


(a)  -  3e(a)  -  3gag  ‘'(a)(a“g)  3gov  "(a)a 


V 


V 


6L_  ,  *  3  L_ 

~2  (a)y  „  ~2 

e  =2 

n  a 

yv 


6e 


(a) 


98 


v 


gfe£-»l  _  .  a*  1 


■ 


31. 


which  is  symmetric  in  y,V  .  The  spin  equations  are  therefore 


(a)[U  6lil  _  n 
e  — 7 — r-  -  u 

x  (a) 


v] 


6l. 


(a)  ~  (a) 


9l  9l 

~1  -  (-^t) 


9L 


oe  i  -v  a e  oe„ 

V  X  V  V 


V 


V 


V 

3ev  3ev  3ev  3ev 


1.  e. 


6L.  9L 


v 


9e 


V 


*  9L 

1-  ( _ L. 

6s  ~*(a) 

3ev 


)) 


3^1  , 1 —  X.  3i^g 

->)  (^~8  n  u  }  |  X  +  .  (a)  n  L: 
de  9e 

V  V 


The  last  term  gives  a  term  symmetric  in  y  and  V  and  disappearing  on 
antisymmetrisation.  So,  assuming  that  the  number  of  particles  is  con¬ 
served,  i.e.  that  Ny  1  =  (/-g  n  u^) 1  =  0  ,  then  it  follows  that  (3.2.11) 

~  I  y  I  y 

and  (2.3.6)  are  equivalent. 


Since  the  Einstein  tensor  satisfies  the  contracted  Bianchi 
identity  GyV|^  =  G  >  the  total  energy  tensor  must  also  satisfy 
T^V|^  =0  if  the  gravitational  field  equations  (3.3.15)  are  to  be  consis¬ 
tent.  One  may  simplify  equations  Ty^|^  =  0  using  the  spin  equations  and 
non-gravitational  field  equations  to  obtain  equations  of  motion  for  Py 
(cf.  the  derivation  of  the  geodesic  equation  for  "pure"  dust  from 

TyV,  =0  with  TyV  =  p  uy  uV)  . 

|V 

However,  a  much  easier  approach  is  to  "vary  the  world  lines"  in 

the  action  integral,  the  problem  reducing  to  a  variation  of  a  one  parameter 

_  „y 

action  integral  of  the  single  particle  type.  The  integral  curves  of  N  , 


. 

' 

' 


' 
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given  by  equations  3x^ (a1, A) / 3 A  =  = 

world  lines  x^  =  x^*(a1,X)  where  a1  , 

the  world  line  (u^  3  a1  =  0)  and  A 

1 

the  curves. 


u^  ds/dA  ,  define  the  particle 
i  =  1,2,3  ,  are  constant  along 
is  an  arbitrary  parameter  along 


/-g  d  x  =  V^g 


9x  '  d3a  dA 


3 (A, a) 


^  eyaBY  vW  E(l)  E (2)  E(3) 


f'L,  d3a  dA 


(3.4.1) 


(E^L^  =  3xa/3a1  ,  v^  =  3x^/3A  ,  is  the  Levi-Civita  permutation 

a 

symbol.)  The  are  tangential  to  the  A  =  constant  surface.  If  n^ 

/  3 

is  the  unit  normal  to  this  surface  and  /  g  is  the  square  root  of  the 
surface  three-metric  then 


a 


8  VgY  E(l)  E(2)  E(3)  *  8  % 


/ -  A^  V1 

v-g  d  x  =  -n  v 


X3g  d3a  dA  =  -n  v^  dj  dA  (3.4.2) 


r  /  3  3 

where  d^  =  /  g  d  a  is  the  3-area  of  section  A  =  constant  of  the 

i  3  v 

infinitesimal  tube  (a  , d  a)  .  The  numerical  flux  across  d^  , 

-n  u^  nl:  dj  ,  is  constant  along  the  tube  since  (n  u^)|y  =  0  •  Denoting 


1! 


3  i  3 

this,  the  number  of  particles  in  the  tube  d  a  ,  by  N(a  )  d  a 


N(ai)d3a  =  -n  u^  n  d£  , 


y 


y 


(3.4.3) 


then  by  (3.4.2),  (3.4.3),  and  the  homogeniety  of  L  , 


L. (Na, . . . )d^x  =  L  (ua,...)  n/-g  d^x 


■ 


■ 


■ 
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L  (ua, . . . )n(-n  vM)d£  dA  =  L- (vU, . . . )n(-n  uM)  dj  dA 

-L  y  j-  y 


=  L^(va,...)  N  d3a  dA  . 


L-,  d  x  = 
~1 


N  (a"S  d3a 


L^(v^,...)  dA 


(3.4.4) 


Keeping  the  applied  fields  \pA  ,  g  fixed,  we  extremise  the  action  inte- 

A  yv 


y 


gral  I  for  arbitrary  variations  in  the  flux  field  N  (a  variation  in 

inducing  a  corresponding  variation  in  the  integral  curves  of  ) 

(a) 

subject  to  6e  =0  and  dN  =  0  .  The  condition  that  N  is  constant 

J  a 

3  3 

simply  requires  that  the  number  of  particles  N  d  a  in  a  tube  d  a 

remains  fixed  on  variation.  Since  L0  and  R  are  not  varied  we  have 

~2 

by  (3.4.4) 


dl  = 


N(a1)d3a 


d  L1(vy,. . .)  dA 


(3.4.5) 


and  resultant  equations  are  exactly  as  in  Chapter  Two,  §2.3. 


§3.5  Special  Case;  Maxwell  Einstein  Field. 

Specialising  the  results  of  this  chapter  to  the  media  in  inter¬ 
action  with  an  electromagnetic  field,  we  take  the  fields  to  t^rie 

four  vector  potential  . 

Making  the  same  assumptions  for  the  dependence  of  L  upon 

as  in  §2.4  for  the  particle  Lagrangian,  we  take  L  to  depend  on  A^  only 

06 

through  a  term  e  A^  N  and  to  depend  upon  only  via  =  A^|^  “ 

Ae|a  • 


■ 

.  '-.a 


«.’•  3'.  J  .  3 
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We  have  at  once  the  electromagnetic  field  equations 


0  = 


6L  3L 


3l 


5A  3  A  oA  i  d  3 

a  a  a  3 


i.  e. 


(3.5.1) 


(X  ^  m 

defining  current  J  and  antisymmetric  tensor  H  by 


a  3 


1 -  i 

/-g  J 


a 


9L 

r*J 

Ja~ 

a 


=  e 


a  / —  < 

N  =  v-g  e  n  u 


a 


(3.5.2) 


/- g  H°^  =  -8tt 


3L  3L 

=  4tt 


3F  Q  '™  3A  \Q 

a3  a  3 


(3.5.3) 


a  ft 

Antisymmetry  of  H  and  equations  (3.5.1)  imply  the  conservation  equa¬ 


tion  J01!  =0  for  charge,  and  J1*  =  e  n  u*  ,  J  i  -  0  and  (n  u  )  i  -  0 


„a 


a  ,a 


a. 


a 


a 


yv 

together  imply  de/ds  =  0  .  The  canonical  energy  tensor  t^ 
field,  given  by  (3.3.13),  is 


a 


for  the 


yv  ,  .-1/2  _  yv  „  .-l  .  y  uotv 

c»)  =  (-g)  ~2  g  -  (4w)  V  H 


(3.5.4) 


and  the  field  spin  flux  (see  definition  at  top  of  page  28)  is 


(4*rl  H“V  (-sP[T6a])  Ae 


=y^-g  (2ti)_1  HV[yATl  =  v^i  (2tt)— 1  A[V]V 


(3.5.5) 


This  implies 


■ 


* 
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,yxv  ,  _vxy 


-1  u  TTXV 


s7;;  +  s;:r  +  =  /-g  (2tt>  x  am  h' 

(ip)  ~(ip)  ~0j;)  6 


(3.5.6) 


(3.5.4)  and  (3.5.6)  together  give,  using  (3.5.1), 


,  1,  „x-l/2  ,cyxv  vxy  vyx  . 

(W  2(  g)  (S((!)  §(W  )|T 


-a"  JV  +  (4tt)_1  f/  H™  +  (-g)-1/2  L2  g"V 


(3.5. 7) 


Now  defining  a  new,  gauge  invariant,  "kinetic"  momentum  p  as  in  §2.4, 


y 


yv  _  y  v 

p,  =  P  -  e  A  .  Matter  tensor  T  .  N=np  u  is  then 
y  y  (mat)  K 


Tf  ,  -  n  P1*  UV  -  Jv 
(mat) 


(3.5.8) 


Define  a  gauge  invariant  electromagnetic  tensor  by 


myv  -  //  n-i  „  y  TTxv  ,  ,  .-1/2  _  yv 

T  (e.m. )  =  (4,T)  Ft  H  +  (-g)  ~2  8 


(3.5.9) 


Equations  (3.5.7),  (3.5.8)  and  (3.5.9)  now  give  the  final  form  of  the 
gravitational  field  equations  (3.3.15)  as 


T^V  N  =  n  p^  uV  and  S^!TV  .  =  n  S^T  uV  are,  respectively,  the  material 
(mat)  (mat) 


fluxes  of  momentum  and  spin. 


' 


. 


■ 
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If  we  choose  for  L0  the  usual  free  field  electromagnetic  Lagran- 

~2 

— 1  j —  n\)  a3 

gian  -(16tt)  /-g  F  v  F  ,  and  define  a  polarisation  tensor  M  by 


Mae  =  3^/aAg,  =  2  3L1/3F| 


a3 


then  (3.5.3)  gives 


4tt  M 


.a  3 


(3.5.11) 


and 


a 


.a  3 


4tt  (J  +  M 


are  the  electromagnetic  field  equations. 


(3.5.12) 


. 


CHAPTER  FOUR 


HIGHER  DERIVATIVE  COUPLING 


§4.1  Generalisation  of  the  Rosenf eld-Belinf ante  Identity  (1.14). 


In  Chapter  Three  we  considered  a  Lagrangian  L  depending  only 

on  variables  and  first  covariant  derivatives  i  .  We  now  discuss 

A  -  A  a 

rv/  ' 

the  more  general  situation  where  L  is  allowed  to  depend  on  higher  deri¬ 


vatives  of  $ 


A  ‘ 


Since  any  derivative  i  can  be  rewritten,  making  use 

J  A  a. . . .a 

~ 1  1  n 

of  the  Ricci  commutation  relations  (1.8),  in  terms  of  the  symmetrised 


derivatives  $A  ,  *  ,  , . .  .*A|  (a  )  .  >  • 

~  ~ 1  1  n 


R 


a 


8y6 I (a- . . .a  ) 
1  1  n 


we  can  assume  that  our  Lagrangian  is  dependent  only  on  symmetrised  deriva¬ 
tives,  L  =  L($A|a(n))  ,  n  =  0,1,2,...  ,  a(0)  denoting  the  empty  set  and 


a(n)  ,  n  _>  1  ,  denoting  a  symmetrised  set  of  indices  a(n)  =  (a^...oO  . 

A  repeated  set  of  indices  a(n)  will,  in  the  usual  way,  imply  summation 
over  the  respective  indices.  The  use  of  symmetrised  sets  enables  one  to 
deal  with  indices  without  any  worry  about  their  ordering. 

The  first  step  is  to  generalise  the  identity  (1.14),  the  condi¬ 
tion  that  L  be  a  scalar  density.  We  consider  a  scalar  density 


~  ~(<I>A|a(n) 


)  ,  n  =  0,1,2, .. . 


(4.1.1) 


and  make  the  following  definitions: 


37  - 
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a  Aa  .  •  .ot  .  ✓  \ 

LA  E  3L/90a  ,  L  1  n  =  L^(n)  =  3L/3^»a  .  .  „  . 

A  ~  ~  ~  A  (a_ . .  .a  ) 

1  1  n 


(V<  <v 


(4.1.2) 


Aa, . . .a  .  /  N 
,  1  n  Aa(n) 

L  =  L  ~ 

^4JL  * 

/>  rsJ 


y  (-i>m  LA~(n)~(m) 

L»  oj 


m=0 


lg(m> 


(4.1.3) 


UT°  E  (I  B)  a  l  Lf2(n) 
~  p  A  P  n=0 


0 


B I  a(n) 


(4.1.4) 


w 

a  _  ra  T  r  /  ,,n  *  T Aaa(n) 

t  =  6  L  -  I  (n+1)  0.  i  ,  ,  v  ^  L  ~ 


V  P 


n=0 


A  (pa(n) ) 


(4.1.5) 


The  condition  (1.10)  that  L($A)  be  a  scalar  density  is 


3L 


(I.~)  a  $_  +  6Q  L  =  0 


30.  A  p  B  p 


(4.1.6) 


For  a  Lagrangian  (4.1.1),  0^  =  {O^i  ,  n  -  0,1,2,...}  ,  this  is 


OO 


I  LA2(n)  (I 

n=0 


Aa(n) 


B£(n))  °  ,  +  6°  L  =  0  .  (4.1.7) 

' p  B  |  ^  (n)  p  ~ 


Repeated  use  of  (1.5)  gives 


LAa(n)  (1 

La  V-1- 


B8(nh  P  .  _  T  Aa(n)  .  B  a 

Aa(n)  ^p  ^B|_g(n)  ~  A  p  B|a(n) 


TAaa(n-l)  _  _  n  9 

*n  ~  $A|(pa(n-l))  ’  ” 


(4.1.8) 


Substitution  of  (4.1.8)  into  (4.1.7)  gives  identity  (4.1.7)  in  the  form 


w 

l  hA~M  U.  >-a  *r.  •  •  +  t.-  -  o  . 


n=0 


A  'p  B|a(n)  ~p 


(4.1.9) 


' 


>1  ■  ■ 


*  0° 


. 


:-Vv  ^.  (?,I).  in 


' 


39. 


From  (4.1.3), 


ATa(n)  Aa (n)  TAa(n) 

L .  ~  i  =  L  ~  -  L  ~ 

*T 


(4.1.10) 


Therefore  the  divergence  of  (4.1.4)  gives 


uT0  i  =  a,B)  0  l  (a 

P  T  A  ' P  ~  ~ 

n=U 


Aa(n)  TAa(n)  4  LAT«(n)  s 

*  ;  B  a(n)  B  ka(n) 


i.e. 


» 


TO 

U  | 

~  p  T 


) 


a 


P 


co 


l 

n=0 


LAot(n) 


$ 


B|a(n) 


) 


a 


P 


(4.1.11) 


Finally,  the  insertion  of  the  above  into  (4.1.9)  gives  the  generalised  form 
of  identity  (1.14): 


U 


TCJ 


P  I  T 


L* 


(iB)  0  *„  +  ,  a  =  0 

A  p  B  ~p 


(4.1.12) 


§4.2  Tetrad  Variation. 


The  action  integral  is  now 


I  = 


(-(16ir)-1(-g)1/2  R  +  L(fA|a(n))> 


(4.2.1) 


with  =  {Na,e^a\ if.  ,R°k  „}  .  Symmetrised  derivatives  of  any  order  are 
A  ~  a  A  3y<5 

/V  1 

now  allowed  to  appear  in  L  ,  but  we  assume  again  that  no  derivatives  of 

ct  (  ^  ) 

N  occur  and  only  first  derivatives  of  e  are  allowed,  i.e., 

~  a 


are  allowed.  i.e 


' 

•i 

. 

. 


■ 


o 
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^Ala  ^Ala(l)  ^  eot  1 3  ’  ^  A 1 3  ’ R  3y<5|p^ 


(4.2.2) 


and 


^A  |  ot  (n)  ^A|a(n)  ,R  3y5  |ot(n) 


}  ,  n  >  2 


is 


(a) 

For  a  variation  de  in  the  tetrad  field,  the  change  in  L 

Ot 


dL  =  l  L— (n)  d($A  J 

^  U  r 


n=0 


A  a(n) 


(4.2.3) 


From  (1.7)  the  variation  d($  i  ,  for  n  >  1  is  given  by 
N  A  a(n)  — 


d($A|a(n)}  (d(<l>A|a(n-l)))|a  +  (IAa(n-l) 


BB(n-l),  oA  d  rP 


}p  $B|g(n-l)  d  oan 


(where  a  here  is  "symmetrised"  along  with  the  other  o^(n-l)).  (4.1.8) 


n 


then  gives  (for  n  _>  1  ) 


Aa(n)  f  x  Aa(n) 

L —  d($A  |  f  \ )  =  L — 

~  A  a(n)  ~ 


(d($  |  ,  .  J)  | 

A  a(n-l)  a 


n 


+  L5a(n)a  S)  a  a.  ,  ,  s  d  rp 

~  ^  A  'p  n  a(n-l)  aa 

rst  ~  rsj  \  II 


-(n-l)L^°S(n‘2)an  $ 


A|(pa(n-2))  d  raa 
-■ 1  n 


(4.2.4) 


Substitution  of  the  above  into  (4.2.3)  and  an  "integration  by  parts  gives 


co  A6 ,06  (n) 

dL  =  L~  d$A  +  (div)  -  7  L  i„  d($ 

r+j  />«»  A  U  rs* 

~  n=0 


8  V‘A  a(n) 


) 


a  101)  asvlg  audi 


. 


' 
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/T  Bn  a  v  ATa (n)  A 

+  aA~)p  l  h~  * 

~  ^  n=0 


,  d  rp 

B  a(n)  ax 


rs_/  l 


l  (n+1)  $ 

n=0 


L^a(n)  p 

A|(pa(n))  ~  ax 


(4.2.5) 


A$  a(n) 

Repetition  of  the  above  procedure  (we  now  have  -L  j „  replacing 

I 

L43:(n)  in  (4.2.3),  (4.2.4))  will  give 


$  ,  arp 

a  '  V‘A  'p  **  B|a(n)  ax 


dL  =  (div)  +  L~  d$  +  (IA~)na  I  L*  ~ 


/N/  ^ 


V  /  1 1  \  /K  TAaxa(n)  rp 

)  (n+1)  $.  |  .  .  N .  L~  ~  dr 

Lr.  A  (pa(n) )  ax 

n=0  1 


/v  1  r*o 


(4.2.6) 


Define 


xa  _  v  ,  _ Aaxa(n)  .  to  xa  xa  .  . 

~  p  n  ^  1)$A|(pa(n))  ~  ’  ~*p“~p*~p 

n=U  ~  ^ 


then  comparison  with  (4.1.4)  yields 


dL  =  L~  d$A  +  uTa  dFP  +  (div) 

~  A  p  ax 


(4.2.8) 


A(a^. .  .cO 

Define  the  gravitational  multipole  moments  Q  and  associated 

quantities  as  special  cases  of  general  definitions  (4.1.2),  (4.1.3) 


9A BY6  =  9A  =  3i/3RA  -  3L/3RXgY6 


(4.2.9) 


.  ,  .  Aa, . . .a 

9  2<n)  =9  n  =  %/3RA|a(n) 


(4.2.10) 


9^  =  9^  l  (-Dm9Ag(m) 


m=0 


l£(m) 


(4.2.11) 


4 


^  J 


\ 
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The  set  d$A  is  {de^jdR*^}  ,  giving 


L~  d$A 
A 


6L 


4(a)  y 

'y 


de (a)  +  2  3y6  dR01 


sa* 


3y6 


and  therefore  (4.2.8)  as 


dL  = 


-P—  de^a^  +  U^Q  drP  +  3  3Y<5  dRa  .  +  (div)  .  (4.2.12) 

(a)  y  p  OT  PY o 


'y 


xa 


3y6 


This  equation  is  identical  to  (3.2.7)  except  that  ^  and  3a* 
replace  UTa  and  Q  3Y^  ,  and  a  calculation  identical  to  that  of  Page 

f  ~  p 

24  gives,  as  the  generalisation  of  (3.2.10),  the  gravitational  and  spin 
equations : 


-(8ir)  1(-g)1/2  GyV 


6~  (a)y  my[VT]  nv[yT]  TJT(yvk 
-Jo)  e  +  (K*  +  ;|t 


yX 


(4.2.13) 


Since  only  one  term  is  not  symmetric  in  y  and  V  ,  as  in  Chapter  Three, 
the  spin  equations  are 


(a)[y 

e 


6l 


0 


Again  the  Belinf ante-Rosenf eld  identity  is  used  to  eliminate 
in  (4.2.13).  For 


(4.2.14) 


6l 


6e 


(a) 


V 


L 


L($ 


A|a(n)  ^ 


ro 


(V» 


(4.2.15) 


' 

■ 


■ 


/j.rqfc  bru 
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\|a(n)  {^A|a(n)’R  6y6|a(n)}  »  nl2  > 


(4.2.16) 


/VI  ■  /VI 


we  have  (see  (4.3.3)) 


iA  (I  2, 0  ♦  =  J~  N° .  J~  Na  6°  - 
A  ;p  £  3np  ~  3Na  ~  p 


6l 


(a) 


(a)  'p 

> 

'O 


+  **»)(IaV  *b  +  2X\0rb 


(4.2.17) 


where  denotes  the  special  case  of  definition  (4.1.3)  with  ^ 

set  equal  to  ip^  .  i.e. 


£*0J0  "  (_1)  (8£/^A|3(m))|3(m)  * 

m=U 


(4.2.18) 


The  non-gravitational  field  equations 


£*(iJ0  0 


(4.2.19) 


and  (4.2.17)  give  the  Rosenfeld-Belinfante  identity: 


UTa  |  +  L~  (I  ~)  a  $  +  t  °  =  0 

~  p  T  v  A  'p  B  ~p 

1  *  rv/  rw 


in  the  form 


A_  e(«)  .  N°  _  !L  N“  S°  +  uTO  +  t  0  +  gj(i  ! VV  (4.2.20) 

6e(a)  P  »NP  ~  3n“  ~  P  ~  plT  ~P  A  P  ^ 

o  ~ 


Substitution  of  (4.2.20)  into  (4.2.13)  gives  the  gravitational  equations 


as 


■ 


r  ;■ i  ii  1- 


•  • 


■ 
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(8TT)_1<-g)1/2  o'”  -  4  SV  gPM  -  4  N“  gWV  +  tyV 


3Nf 


3N 


a 


U[VT]  v[yx]  x[vy]  +ytvu 

*  '  t  ~  T 


^2l(lJV)Y|YX  +  st«AVvR 


B 


(4.2.21) 


§4.3  Lagrangian  Decomposition. 


We  continue  in  similar  fashion  to  the  derivation  of  §3.3, 
writing  L  as  a  sum  L  =  of  a  free  field  Lagrangian  for 

and  an  L^(N^,...)  homogeneous  in  ,  L^(N^,...)  =  n/-g  L^(u^,...)  . 


~  ~l®y>ea  )>ea  ),'|,A|ci(n)’RA|a(n)')  +  ~2  %  I  a  (n)  ^B5  ' 


(4.2.9),  (4.2.10)  and  (4.2.11)  become 


^  -  SA  -  nyCg  3Ll/3RXM  =  n^i  =  n/^  QA 


pAa(n)  _  n/_  3l  /3R  =  n/.g  QA2(rl)  =  n/_g  q  1  n 

~  1  A  a(nj 


qA  -  Qa*M  -  Ai  I  (“l)m(n  QA£(m)) 


m=0 


3(m> 


(4.3.1) 

(4.3.2) 

(4.3.3) 

(4.3.4) 


Aa1 . . .a 

where  Q  n  are  the  gravitational  multipole  moments  per  particle. 

We  continue  by  examining  the  definitions  of  the  other  terms  appearing  in 
(4.2.21),  keeping  in  mind  (4.2.15,  16)  which  imply: 


-  j ;  .Si 


+  . .  *  J  /iuia  £  e a,;  J[  ^nJt3l7w 


' 

0  3-  n  „  .  5x{>\ ,  .1*  PSfii  * 


■ 
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3l  Bl  Bl  . 

T  A  r~  ~  ~  -  T  A 

—  L  «  /  \  %  *\.i.  “i-i 


BL  Bl 


Br, 


’  ^a)  '  ^ ' =<w  ’  ^a6Yfi  "A 


=  qa> 


Aa 

L~  1  =  (- 


Bl 


Be 


(a)  5  dip 


Bl  Aa,  Aa 

=  L,.l  ,  Q  b 


a  a. 


A|a1 


(ip) 


I  \U  )  r^j 


Aa, . . .a  Aa, . . .a 

T~  1  n  _  r  1  n  _ 

~  (ip)  “  Bip 


Bl 


A|a(n) 


,  qA2(n)}  (n>2) 


l~  =  { 


BL  6l 

r-*j 


.  lL„,  .  gt> 


*  V  ’  6e(a)  ’  ~*(W  ’  ~* 
~  a 


Aa  Bl 


Aa.  Aa. 


T~  J-  _  { _ _ _  t  o  1> 

k*  ~  3  (a)  ’  ~*(\p)  ’ 

dea|a 

it2(n)  ■  <£$?  •  5*S(n)}  <*» 


Definition  (4.2.7)  and  (4.3.10)  gives 


TO  V  /  -i  \  i  t  ■ 

~  P  =  n=0  (  ^Al  (pa(n) )-*(*) 


Aoxa(n)  +  j-  (n+1)R 


,Aaxa  (n) 


n 


=0 


A  |  (pa(n)P* 


=  +  yTa 

~(\p)p  ~(R)p 


From  (4.3.9),  (4.3.10)  we  obtain 


U 


xcr  _  B.  a 


(  ,  u  j  jAmfa)  # 

~  M  n=0 


TO  TO  .  ,  r3  ~~  (a) 

=  ~0|>)p  +  ~(R)p  (_6P  6o)  ^(a)  e3 

'ax 


B  |  a  (n) 
Bl 

r> 

Be' 


U 


to  _  B.  a  r  Axa (n) 

(«P  "  (  A  }p  *B|a(n)  ’ 


(4.3.5) 

(4.3.6) 

(4.3.7) 

(4.3.8) 

(4.3.9) 

(4.3.10) 

(4.3.11) 

(4.3.12) 


• 

. 
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TG 


~(R)P  (IA  2* 


Axa(n) 


^IgCn)  - 


(4.3.13) 


Since 


3L 


(a) 


(a)  p 
'a  |  t 


8Li  /  \ 

~1  _ (a) 

,  (a)  p 
>e  i 
a  x 


9l 


=  -nv-g 


1  e(a>  uT 


(a)  'p 


(4.3.12)  gives 


uT[opI  =  i^[pp]  +  u^?p]  -  i  nt  sp0 

~  ~  (ip)  ~(R)  2  ~ 


(4.3.14) 


where  S^°  is  the  spin  angular  momentum  per  particle. 


Defining  the  spin 


flux 


SPTV  E  2  yVtyx] 


(4.3.15) 


then  (4.3.14),  (4.2.7)  give 


Sf V  =  2  U 

7v 


•V [ y  x ]  +  2  yVtyx]  _  2  yv[yx]  +  gyx  Nv 


(ip) 


(R) 


(4.3.16) 


Defining 


SWTV  =  2  uVtpt]  +  2  dv[wt]  _  2  yvIpt] 
~  (ipR)  ~(ip)  ~(R)  ~ 


(4.3.17) 


gives  the  spin  flux  in  the  form 


■  ^VR)  +  ^  sv 


(4.3.18) 


S^T  NV  and  a  field  part 


a  sum  of  a  matter  part 


which  contains 


W” 


"t 
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contributions  from  the  curvature. 


Since  L  is  a  function  of  e 


(a) 


(a) 

a|  3 


and  L,  is  a  function  of 
~1 


e^  ,  a  calculation  identical  to  (3.3.8,  9,  10,  11)  of  Chapter  Three 

LX 


will  give 


3L  „ 

~  a 
N  =  L 


a 


9l 

~  py 

—  N  g 


3N 


u  v 
P^N  + 


1 

3L 


l_e(a)|UNV 


3n 


* 


(4.3.19) 


(4.3.20) 


a 


where  is  the  canonical  momentum  per  particle.  The  final  term  of 

(4.2.21)  to  examine  is  t^  +  Y  .  Using  (4.1.10),  (4.3.9,  10) 


.TV 


+  t  v  =  6V  L  -  l  $ 

I  1  I  /v  Xj 


y  t  ~y  y 


n=0 


Ava  (n) 

A|  (ya(n))  ~ 


(4.3.21) 


(L1+L9)6 

~1  ~2  y 


v  (a) 


9L 


-  e 


a|y  3e(a) 

a  v 


n=0 


V  i  Ava(n)  y 

^  \|(ya(n))  ~*($)  "  A |  (ya(n) ) 


Ava(n) 


n=0 


(4.3.22) 


with 


-  e 


,  \  3L 
(a)  _ 2^ 

o|y  3  (a) 

a  v 


.  e<«)  _?h  Nv 

ay  ^:(a) 


9e 


a 


(4.3.23) 


Define  a  canonical  energy  tensor  for  the  fields  ip.  by  ([14]) 


v 


^  SKtJ/)  ~2  6y  “  J;0  ^a|  (ya(n))~*(ip) 


Ava(n) 


(4.3.24) 


(4.3.19,  20,  22,  23  and  24)  then  yield 


. 


•  '  \ 


v  ;  .  >J  it  1C  10  i-  1  M  ft  .  -it  •'  e  .:j13<h1 


-  4  .  /  \  *  H**x 


US  bff*  CS  ,SS  ,0S  , *!.£.*) 


48. 


3L 


9N 


a 


Na  g^V  +  tyV  +  YTV^ 


T 


U  V 
=  PM  N 


+  /- 


g  t 


yv 

010 


oo 


i 

n=0 


qAva(n) 


RA| (ya(n)) 


and,  finally,  substituting  the  above  expression  into  (4.2.21)  gives  the 
final  form  of  the  gravitational  field  equations: 


-W-V™  -  T*”  5  »  PV  +  +  +  +  S?1') 


T 


+  (-g)  1/2  (-  I  Q 


00 


n=0 


Ava(n)  p 

*  Ka| (pa(n) ) 


*  *  s^V  *  si  <■.*>"”  V 


(4.3.25) 


In  summary,  the  above  equations  generalise  the  equations  (3.3.15)  of 
Chapter  Three,  the  total  energy  tensor  is  found  to  contain: 


(i)  A  material  energy  momentum  tensor  ^^at)  =  n  uV  where 


P^* 1  is  the  canonical  momentum  per  particle  given  by  (2.3.5), 


(ii)  A  canonical  energy  tensor  for  the  fields  given  by  (4.3.24), 

(iii)  A  contribution  from  spin,  the  spin  flux  being  the  sum  of 

a  matter  part  ~^at)  =  ^“8  n  S^T  uV  and  a  field  part 
that  includes  contributions  from  curvature. 


(iv) 


Contributions  from  gravitational  multipole  moments. 


. L  ...  .  *J  -  ;  -  io 


. 
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CHAPTER  FIVE 


DIRAC  EQUATION 


§5.1  Discussion. 

The  possibility  of  obtaining  classical  relativistic  equations 
for  spinning  particles  starting  from  the  Dirac  equation  has  been  explored 
in  papers  by  several  authors.  The  main  attempts  have  been  directed  at 
deriving  the  Bargmann,  Michel  and  Telegdi  equations  [7],  by  making  suit¬ 
able  definitions  for  spin,  four  momentum,  etc.,  either  directly  in  terms 
of  spinor  entities,  or  as  expectation  values.  More  recently,  Suttorp  and 
de  Groot  [9]  derived  from  the  Dirac  equation  operator  equations  having 
the  same  form  as  their  classical  equations  for  a  composite  particle. 

The  two  main  approaches  have  been: 

(i)  to  derive  operator  equations  (see  [12],  and  also  [9]), 

(ii)  to  apply  the  W.K.B.  approximation  to  the  Dirac  equation  (or  the 

squared  Dirac  equation)  [13]. 

In  (ii)  the  lowest  order  equations  have  similar  form  to  the  Bargmann, 

Michel  and  Telegdi  equations. 

None  of  the  references  really  seems  to  answer  the  question  of 
how  one  should  define  relativistic  objects,  such  as  tensors  describing  spin 
and  momentum,  in  terms  of  quantum  mechanical  entities.  This  same  diffi¬ 
culty  arises  in  §5.3  where  exact  equations  are  quite  easily  obtained  from 
the  general  relativistic  Dirac  equation,  equations  having  the  same  formal 
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appearance  as  equations  (2.4.1)  and  (2.4.2).  Section  §5.2  will  briefly 
describe  the  necessary  notation  and  four-spinor  calculus. 


§5.2  Four-Spinor  Calculus. 

Let  Y  (a  =  1.....4)  denote  a  fixed  set  of  Dirac  matrices 
a 

satisfying 

Y(aYb)  =  nab  =  d±ag  t-1*-1’-1’1)  *  (5.2.1) 

If  Y+  denotes  the  hermitian  conjugate  of  Y  >  then 

Y+  =  “Ya  \  =>  Y4YaY4  =  Ya  (a=1>-*-’4)  (5.2.2) 


Define 


5  Y[aYb] 


(5.2.3) 


then 


Y  a,  -  a.  y  =  ^  n  r  J  n 

'a  be  be  a  a[b  c] 


(5.2.4) 


If  ooak  =  -  a)ka  is  an  arbitrary  set  of  scalars  (lxl  matrices)  then 


(5.2.4)  implies 


cl,  N  be 

co  y  =  t  (Y  -  a,  y  )  w 
a  ’c  4  'a  be  bc'a 


(5.2.5) 


(^0 

For  an  o.n  field  of  tetrads  e  (x)  : 


' 


■ 


..  Y" 


■ 
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n 


ab 


■  Wx) 


(5.2.6) 


the  Ricci  rotation  coefficients  are 


a  (a)  a  3 

Y  be  ~  ea| 3  6(b)  e(c) 


Note 


Yabc  Y[ab]c 


(a)  a 

ea|3=-Y  be 


(5.2.7) 


(5.2.8) 


(5.2.9) 


Field  of  y  (x) 
-  1  y v 


For  fixed  numerical  matrices  v  define 
-  i  a 


Y  (x) 
'y 


(5.2.10) 


(5.2.1  and  6)  imply 


Y(yYv)  8yv 


(5.2.11) 


(5.2.9  and  10)  imply 


(a) 

Y  i  E  Y  e  I 
'y  V  a  y  v 


a  (b)  (c) 

-  -Ya  Y  be  ey  ev 


Write  this  as 


■ 

. 


*• 


■ 
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then 


e,.*^  y  I  dxV  =  7-  (y,  a  ,  -  a  y  )  go 
(b)  'y  V  4  b  cd  cd  b 


cd 


by  (5.2.5) 


i.e. 


y  1  =  -  (y  r  -  r  y  ) 

'y  v  y  v  v'y 


(5.2.12) 


_  _  1  cd  (f) 

rv  -  ~  4  acd  Y  f  v 


(5.2.13) 


Define 


V  y  =  y  1  +  [y  ,r  ] 

v'y  y  v  y  v 


(5.2.14) 


then 


V  y  =  0 
v'y 


(5.2.15) 


The  covariant  derivative  of  a  spinor  ip  is  defined  so  that 

y^V  gJj  +  =  0  is  invariant  under  co-ordinate  transformations,  and  also 

invariant  under  arbitrary  tetrad  rotations 


ei$a)  =  A\  (x)  eyb) 


(5.2.16) 


which  can  vary  from  point  to  point.  From  (5.2.16)  and  (5.2.10) 


Yy  (x)  =  Ay  (x)  Yv(x) 


(5.2.17) 


where  =  eV(a)  A\  e<b)  .  Since  y^Y^  gyv 

there  exists  a  4x4  matrix  S(A)  such  that 


by  Pauli's  theorem, 


'  . 


■ 


. 
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Yp  =  S  Ym  S’1  .  (5.2.18) 

Postulate  that  under  tetrad  rotation  (5.2.16)  ^  transforms  as 

(x)  =  S(x)  ip(x)  .  (5.2.19) 

(This  ensures  that  for  constant-tetrad  rotations  (e.g.  special  relativity) 

3  ^  =  0  is  invariant . )  Then 

Yy,3  r  =  S  yy  s"1  3y(S^) 

=  S(Yy3y^  +  Yy(S_13yS)^) 

=  S  Y^O,^  +  (S*13  S)ij;) 

M  r* 

Now,  the  transformation  law  for  T  under  (5.2.16)  is 

r1  =  s  r  s"1  +  (3  s)s_1  .  (5.2. 20) 

y  y  y 

This  can  be  computed  from  (5.2.16,  13  and  7),  or,  more  easily,  from 
Y  I  =  -[Y  ]  and  (5.2.18).  Hence,  defining 

'y  |  v  L  'y  v 

V  li )  =  3  -  T  ip  (Tp=4xl  spinor)  (5.2.21) 

yr  yr  y 

we  achieve  invariance  of  Y^  ^  =  0  .  For  the  Pauli  conjugate  spinor 

1jj  E  \p+y  ,  define  V  acting  on  Pauli  conjugate  ip  by  V^OJj)  =  (vy^)  • 


V  \p  =  3  ^ 
\iy  y 


*  ru  Y/ 


Then  (5.2.21)  gives 


:  >,v 


' 


. 

.  v.  ■;  -  ■}  £  =  *  7 
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=  V  -  *  (VyV 


■  V  +  *  ru 

by  (5.2.2,  3  and  13).  Hence,  for  lx4  spinors 


V  *  =  3^(J)  +  (p  ((j>  =  1x4  spinor)  .  (5.2.22) 

For  tensor-spinors  (4X4)  such  as  y^1  ,  etc., 

v  TKX  =  TKA|  +  [TKA,r  ]  .  (5.2.23) 

y  |y  y 

All  these  covariant  derivatives  transform  as  tensors  under  co-ordinate 
transformations,  and  are  covariant  under  tetrad  rotations  in  the  sense 
(V  ip)1  =  S(V  ip)  ,  (V  4>)f  =  (V  <j>)  S'1  etc.  The  rule  (5.2.23)  is  dictated 

y  y  y  y 

by  the  requirements 


(i)  that  the  product  rule  hold  for  V  ; 

r4, 

(ii)  any  spin  tensor  TkA  =  ip  (p  TKA  or  sum  of  such  terms; 

(4X4)  (4xl) (1x4)  (lxl) 

(iii)  for  plain  tensor  TkA  ,  V^TkA  =  TKA 1  . 


Identity  (5.2.15)  is  unaltered  by  replacing  r  ,  by  , 


where  +  i  I  .  (Factor  i  is  included  to  preserve  identity 


y 

(1X1) 


Y  r  Y  =  -r  (A  real  vector).)  The  Ricci  commutation  relations  take 
T4  y  4  y  y 


y 

the  form 


(V  V  -  V  V  )ii  =  $  ip 
v  v  y  y  v'v  yv  y 


(5.2.24) 


. 


Sfii  .  (  >  t<; .  -TO ft.*  ,3  loH 


..  fiua-al  j  'l'  Vd 


•4  li  Ic  <  sUjS  bo?q  »rl3  Ibrtl 

■ 


. 


where 


(5.2.25) 


3  r  r  - 
[y  v] 


[r  ,r  ] 

y  v 


j  R  R  aa3 
4  yva$ 


+  2i 


3  r  A  . 

[y  v] 


(5.2.26) 


§5.3  Dirac  Equation  in  General  Relativity. 


The  general  relativistic  Dirac  equation,  including  an  anomalous 
magnetic  moment  term  is 


i  yyV^  +  (m/h)ip  +  i  g  FUVW  ip  =  0  .  (5.3.1) 


The  hermitian  conjugate  equation  is 


-i(V)J«V+  +  (m/h)l|)+  -  i  g  yVV+  -  0  . 

Multiplication  on  the  right  by  -Y4  and  insertion  of  Y4Y4  =  1  in 
appropriate  places  gives 


+  i(V^)+Y4Y4Yy  y4  ~  (ra/h)ip  +  i  g  Fyv^Y4YV 


y 


y4y4y  y4 


=  0 


and  use  of  (5.2.1)  and  (5.2.2)  simplifies  this  to  the  Pauli  conjugate  of 
(5.3.1) 


i(V  Wy^  -  (m/h)i|i  tigF^t  yV  =  0  .  (5.3.2) 


Define 


. 


■ 
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T^V  is  a  tensor  (1*1  matrix)  so  V  T^V  =  T*JV .  .  Using  (5.3.1),  (5.3.2) 

M  j  p 

and  (5.2,24,  25,  26)  gives 


Vv  Twv  -  Vv»vV« 

=  (yovV*  +  ?(vvvV« 

=  (Vv«YV(Vy4.)  +  ^(VyVv  +  $yv)YV  <P 

=  ((m/ih)*  -  g  Fag  *Y6y“)  (V^) 

+  l(iVy(-(m/ih)l|;  -  g  I'agYaY6i(') 

+  *<-  1  ^vae0^  +  1  • 


Cancellation  of  terms  gives 


V 

V 


\  Ryvag0J«aV«  +  Fyv»  i  Yv«  -  g  Fag||J(*aa%)  .  (5.3.4) 


Equation  (5.3.4)  has  the  same  form  as  (2.4.2).  Defining 

euV  =  i  \p  y^ip  ;  M01^  =  -2g  ip  o^^ip 


and 


ca3v  _  1  t  a3  v. 

S  =  j  \p  o  y  ^ 


gives  the  R.H.S.  of  (5.3.4)  identical  to  the  R.H.S.  of  (2.4.2).  For  spin 

ry  O  >  \ 

equations,  we  take  the  divergence  of  S  p 


■ 
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VA  s°^A  =  ~  V^a^yNp) 

=  \  (V^a^y'Sj/  +  tyoa^y\xi) 

=  ^  (Vx^)  (YXaa^ -4gA[ay^b^  +  y  ^  aa^yAV^ 

1  ,  m  -7  ^  T  v  yN  a3, 

-  2  (ih  *  - g  V  ^ Y  )a  ^ 

-  2(V[a»YS1lH-  aae(-  ^  ip  -  g  Fyy  yW) 

—  a6 

Cancelling  the  ip  O  ip  terms,  and  noticing  that 


v  y  a3  .  a3  y  v  D_  [a  3]y 
F  Y  y  a  +  F  a  y  y  =  8  F  a 
yv  1  yv  '  '  y 


gives 


V  sa$A  =  _2(VC<V)y3]H»-4  g  F  [°V/]lV  . 

A  y 


Defining  Pa^  by  p0^  =  -  (V^y^  gives  the  above  equations  as 


1  Vx  sa$X  =  _  F[a  M33y 


(5.3.5) 


which  may  be  compared  with  spin  equation  (2.4.1).  The  question  of  inter¬ 
pretation  of  (5.3.4)  and  (5.3.5)  will  not  be  pursued  further,  except  for 

remarking  that  a  quantum  theorist  might  demand  definitions  of  uA  ,  MyV 

yv  .  / 

etc.,  as  expectation  values  rather  than  1X1  spinors.  T  is  a  (canon¬ 
ical)  energy  tensor,  and  one  may  notice  that,  in  geodesic  co-ordinates, 
integration  of  (5.3.4)  over  3-space  (i.e.  taking  expectation  values  of 
the  operators  aa^yV  ,  i  yV  ,  a°^)  gives  for  the  LHS  of  (5.3.4) 


' 


. 
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ThV 

I  i  d  x  = 
V 


TyV  d3x 
»V 


T^1  .  d3x  + 

9 


.  d3x  (i  =  1,2,3) 
>4 


(9Ty4/9t)  d3x  (by  Gauss's  Theorem) 


9P 


y 


9t 


(Py  s 


Ty4  d3x) 


This  is  valid  in  special  relativity  and  may  be  assumed  in  general  relativity 
3 

for  small  d  x  . 


■ 


' 


APPENDIX  ONE 


CO-ORDINATE  TRANSFORMATIONS 


Consider  a  general  co-ordinate  transformation  x13  -  x^(x^) 


Define  =  8x^/8xa  .  Since 


9x°  8xs  _  a 

3xS  3*P  =  P  ’ 


the  inverse  matrix  is  X  =  Bx^/Sx13  . 

P 


A  relative  tensor  $^(x)  will  transform  according  to 


Vx)  =  \B  (XV  Vx)  • 


(A. 1.1) 


Here  A  denotes  the  collection  of  indices,  contravariant  and  covariant, 
of  <f>^  ,  and  repeated  capital  indices  indicates  summation  over  each 
index  of  A  .  All  indices  will  run  from  one  to  four,  unless  otherwise 
indicated.  For  a  scalar 


g 

$A  =  (f>  ,  we  have  A^  =  1  , 


for  a  contravariant  vector 


$A  =  <f> a  we  have  =  X°^  , 


for  a  covariant  vector 


B  -lb 

$  =  <J)  a  we  have  A^  =  X 


and  generally  for  a  relative  tensor  of  weight  w 
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CL  .  ..CL 

$  =  $  1  m  (x)  , 

A  a  . . .a 
1  n 


A 


B 


A  is  given  by 


B  a,  a  _b-  .b 

A  B  =  xRx  ...  xRm  X  1  1  ...  X  1  n  | XP  |  w 
A  3,  3  a..  a  1  a1 

lml  n 


(A. 1.2) 


B  O 

Define  infinitesimal  generators  (1^  by 


a 

'  A  'p  p 

dX 


B 


(evaluated  at  X*"*  -  613  ) 

0  0 


(A. 1.3) 


B  CT 

The  (I.  )  '  s  may  be  constructed  explicitly  from  recursion  formulae: 

A  p 


For  a  relative  invariant  (p  of  weight  w  , 


_  ~  a  _ 

i.e.  4>(x)  =  |~ ~|W  ^(x)  ,  we  have  (p  =  | X |  W(£  =  Acp  , 

9xP 


E  rf  (xPa -  O 


3X 


-w| x| -w_1  iliL  cP  =  SP  ) 

11  0  O  O 


a 


dxr 


O 


=  -w  6 


(A. 1.4) 


Differentiating  X0^  X  =  6^  w.r.  to  XP^  gives 


oOt  *0  -13  ,  va  9  .  -13  v  n 

6P  6B  X  Y  +  x  B  (X  Y>  “  0  • 


a 


-1A  9 

Multiplying  by  X  gives  — 

01  ax' 


,  -1A  x  -lA  -la 

P  (x  Y)=-X  PX  Y 


Hence 


— i-  (X  lX  )  ,vp  .p.  =  -  <5X  «a 
3xpa  Y  (X  o~  V  p  Y 


(A. 1.5) 


$ 
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_  _  g 

For  a  relative  tensor  $  transforming  as  <f>  (x)  =  A  $  (x)  ,  a  rela- 

A  A  A  g 

tive  tensor  $  with  an  additional  covariant  index  OL  ,  will  transform 
Act 

as 


*Aa<x>  =  AAaB6<>  Vx)  '  V  X'1Sa  Vx) 


(T  »e.  a  _  __3_  A  B  -IB 
(  Aa  >p  3xp  A  01  (x  a  =  6  A 


i.e.,  (I.  BB)  °  =  (I,B)  °  «B  -  6®  6° 

Aa  'p  A  p  a  Apa 


(A. 1.6(a)) 


Similarly,  for  an  extra  contra-variant  index. 


-r-a .  B  a  ,3,  \ 
*A«  "  AA  X  B  VX> 


.  /TaBN  a  9  ,,  B  N 

••  <Wp  3xP  (AA  x  6  (XPa  =  5Pa) 


.  aBN  a  B.  a  ra  ,  XB  xa  xa 

i.e.  (I.0)  =  (I.  )  oa  +  o.  6 

A3  p  A  p  3  A  p  3 


(A. 1.6(b)) 


rj-  p 

From  the  group  property  of  transformations,  (Y  XM^) 

A  ^(Y)  A  ^(X)  ,  by  differentiation  we  have 

A 


/ T  c  y  B  3  ,r  c  3,T  B  y  B  3  M  fT  B  y  3 
(IA  h  (IC  >a  -  (IA  >o  (IC  h  -  (IA  h  Sa  -  (IA  >a  5X  ■ 


(A. 1.7) 


B  cf 

The  main  usefulness  of  (I.  )  is  that  covariant  differentia- 

A  p 

tions,  and  variations  on  co-ordinate  transformation,  of  a  relative  tensor 


can  be  succinctly  written  down. 
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For  an  infinitesimal  transformation  xP  =  xP  +  £P(x)  ,  X^  =  3xP/3xC  = 
6^  +  ^  Q  .  Therefore  a  relative  tensor  will  transform  according 

to 


Vx)  '  AAB(6Pa  +  O  Vx> 


3A 


B 


-  AAB(0$B«  +  Hr  (xp  -  6P  )  (5P,a>  Vx> 

3x  a  a 


i.e.  tA(x)  =  $A(x)  +  (IAB)p°  (CP>a)  »b(x)  ■ 


(A. 1.8) 


B  O’ 

(I^  )p  can  be  written  down  explicitly:  For  a  relative  tensor 


a.  ...  a 

f  -  4-  1  m  , 

A  a.,  ...  a 
1  n 


„  a  a  b,  - b 

A  B  _  Y  1  Y  m  Y_1  1  Y-1  n  I Y I _w 

A  p.  Pa.  a 

lml  n 


g 

3A  m  a  a.  a  a  a  b 

—5-  =  (  I  Xg1  ...  Xg1  1  (6  V  JXg1  1  ...  Xgm  x  11 
P  .  ,  P,  P.  ,  P  P.  P  P  a 

3X  „  1=1  1  l-l  1  l+l  m  1 


.  X 


,b 
-1  n 


n 


n  ot  a  b  -b.  -  ..b.  .b.  M  -b 

+  £  xj  ...  x„m  x  1  1  ...  x  1  J_1(-x  1  Jx  1  )x  1  J  1  ...  x  1  n)|x| 
j=l  61  Bm  al  J-1  P  aj  aj+l 


a-  a  ,b,  ,b  1  alyl 

+  (Xg1  ...  x  m  x  1  1  ...  X  1  n)(-w)  |x[  W  AJ 
B1  Bm  al  an  9XP 


a 


„  B  ,  V 
A  p  9XP 


m  an  cx .  -  ot .  a.,,  a  b-  b 

r  r.  1  r  1—1  ,  r.  lpCJ  N  r  1+1  r  m  -  1  f  n 

i=l  “l  °i-l  P  °i  “i+1  ^m  1  n 


-w 


■ 
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n  a  a  b, 

-  I  6g  •••  V*. 

j=l  P1  al 


b .  b .  b 

5  3-1  (6  J6  )6  31 
a.  n  Pa.  a 
3-1  J  J+1 


.  6 


n 


n 


(A. 1.9) 


-w  0  .  0  ^ 
A  p 


B  a 

We  can  write  ^|T  conveniently  in  terms  of  (1^ 


a_ . . .a  m  a.  a. ...a.  _$.a. . . . . .a 

*  ,  =  3  $  1  m  +  l  rx  4  1  1_1  1  1+1  m 

A|x  xa...a  ^  3,1  ar.. 


•  •  •  3. 


n 


n  b ,  ot-  •••  •  •  •  oc 

-  y  r  J  4  1  t  m  -  „  4, 

1=1  J  1  3-1  J  3+1  n 


At  A 


.  cl.  -  a.  „  cl.  a  b 

P  ,  V  *1  ^  1-1  ^  1+1  ...  6  m  6  1  ...  6  n 


m  a 

3  4  +  rp  (7  6/  ...  6  x'x(6  x&'j0  )6 

t  A  ax  8X  3±_1  P  6±  31+1 


3_  a 

m  1 


b 

j  1 
a 


n 


n  0 t  a  b  b .  b  b  b  3,  .  .  •  3_ 

l  6.1  ...  5„m  6  1  ...  6  3-1  (6  J6a  )5  J+1  ...  6  n)4  1  m 

.  3,  3a_  a .  -  p  a ,  a  b-  ...  b 

3=1  1  m  1  3-1  j  aj+1  n  1  n 


+  rp  (-w  sV>  <s> 
ax  p  A  B 


Comparison  with  (A. 1.9)  then  gives 


a  o  +  rp  (i  B)  a  oR 

x  A  ax  A  p  B 


(A. 1.10) 


In  geodesic  co-ordinates,  with 


0  , 


*  .  =  ($  .  >  =  o  o.  -  rp  (i  B)  °  $„) 

A  |  yv  '  A|y  ,v  v  y  A  cry  A  p  B,v 


=330.+ 
v  y  A 


rp  (I  B)  c 

ay,v  A  p 


0 

B 


+  *A- 


■ 


■  .  :  ... 


aV 


>> 
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(r 


ay,v 


-  r 


P 

av,y 


a 

P 


=  -R 


ayv 


a.B)  a 

A  p 


So,  in  any  co-ordinate  system, 


$  |  -  <& ,  =  -RP  (I  ®)  0  $_ 

A|yv  A|vy  ayv  A  p  B 


(A. 1.11) 


APPENDIX  TWO 


CO-ORDINATE  INVARIANCE 


We  write  down  the  identities  resulting  from  the  relative  invar¬ 
iance,  under  arbitrary  co-ordinate  transformations,  of  scalar  densities 
dependent  upon  tensor  arguements. 


Basic  Identity.  Suppose  L($.)  is  a  relative  invariant  of  weight  w  , 


L(y*»  =  lS|Wi(Vx)) 


(A. 2.1) 


9 


x 


—  P  P  rP  /  \ 

Making  an  arbitrary  infinitesimal  transformation  x  =  x  +  t,  (x)  ,  to 
first  order  we  have 


r-O  ra 

=  o  -  E  ; 

rx  cr 

1  1  -  1  -  „  ; 

a  o 

|  9x  |  w 

P  ,P 

3xP  >a 

9xP 

=  1  -  w£ 


a 


,o  9 


giving  (A. 2.1)  as 


L($a(x))  -  L($a(x))  =  -w€  ^  L 


(A. 2. 2) 


Noticing  (A. 1.8),  $A(x)  -  $A(x)  =  (IA~)pa  $B(£P}G)  ,  gives 


9L 


3$  Op°  V5”  a>  =  'w,5p  t(sP.o)  for  arbicrary  ZP<X>  ■ 

A  5 


L  therefore  satisfies  the  basic  identity 


3L 

TTT1  (I.~)  a  +  w6a  L  =  0 
9$.  A  'p  B  p  ~ 


(A. 2. 3) 
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As  a  special  case  of  the  above,  take  to  consist  of  a  set  of  tensors 

(N/ 

and  their  first  covariant  derivatives  i  .  If  L  is  a  scalar 


ha| 


a 


density,  w  =  1  ,  (A. 2. 3)  is  then 


3^  (IAB)p°  +  3*^  V  V  *B|  e  +  5P  i  ■  0  • 


Using  (A. 1.6 (a)), 


9l  _  _  9l  „  9L 

rv  B  Cf  R  (J 

-  (I  )  OH - (I  )  0  i  —  - 

3$a  A  >p  3$A|a  A  >p  B|a  3$A|a 


4>  ,  +  6°  L  =  0 

A I  p  p  ~ 


Rearranging  terms, 


( — 2:1  -  V  — )  (i  B)  °  $  +  ( —  (I  B)  0  $  )  | 

(**A  o3*A|a  A  P  B  WA|a  A  P  B,a 


9L 


9$ 


a|  a 


$A  i  +  6  L  =  0 

A|  p  p  ~ 


l.  e. 


S  (I  B)  a  *  +  ua0  ,  +  t  0  -  0 

6$  A  p  B  ~  p|a  ~p 


(A. 2. 4) 


where  we  define 


a  9L  A  9L  6L  . 

A  _  _ 21  t  Aa  _  ~  _ 21  =  T  A  _  Aa 

~  "  9$A  ;  ~  "  9$A ,  ;  6$A  ~  ~  la  ; 

A  A  a  A  1 


TTa a  _  TAa,T  B.  a  a  _  xa  .  TAa 

u  =  L  (I .  )  ;  t  =  L  6  -  OiL 

~p  ~  A  p  B  ~p  ~p  A  p  ~ 


. 
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APPENDIX  THREE 


DEFINITION  OF  SPIN 


Consider  the  simple  example  in  classical  mechanics  of  a  sphere 
rotating  about  its  centre,  with  kinetic  energy  of  rotation  I  w2  and 
angular  momentum  I  $  (I  =  moment  of  inertia,  g  =  angular  velocity). 

Let  e,  .  ,  i  =  1,2,3  ,  be  an  orthonormal  triad  rigidly  attached  to  the 

~(i) 

sphere  with  origin  at  the  centre.  Since 


^•  =  HXS(1)  (i  =  1,2,3) 


de/>N  de,.. 

-1--  •  -----  =  (w*e  (  v  )  *  (w><e  (  . .  ) 

dt  dt  ~  ~(i)  ~  ~(i) 


, Einstein  summation  ^ 
convention  used  here 


-  '  S(1)  <%) 


3w*w  -  (w’e, .  >.)  (w*e,.N) 


2  2  2  2  2 
3w  -  (w^  +  1*2  +  w3^  =  2w 


The  kinetic  energy  of  rotation  is  therefore 


12  1  T  d~(l)  .  d~(i) 

2  IW  “  4  I_ di - dt 


(A. 3.1) 


Since 


de  . 

x  e/jX  =  e/.Nx(e/.  vXw)  =  e  . . ,  (e, . *  *w)  -  w(e,.x*e,  *) 
dt  ~(i)  ~(i)  ~(i)  ~  ~(i)  ~(x)  ~  - (i)  ~(i) 


=  w .  e .  -  3w  =  w  -  3w  =  -  2w 

A  n*j  -i  r>*>  r*j 
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The  angular  momentum  Iw  is  therefore 


T  =  A  T  v  d~(i) 
1  ~  2  1  ~(i)  dt 


(A. 3. 2) 


(a) 

In  terms  of  an  orthonormal  tetrad  e^  ,  (a, a  =  1,...,4)  with 

Uv 

e^4)(s)  =  (0,0, 0,-1),  (s)  =  (e(i),0)  ,  i  =  1,2,3  and  with  t 

replaced  by  proper  time  s  ,  since 


=  de^/ds  =  (0,0, 0,0)  ,  e^  =  (de^^/dtjO)  i  -  1,2,3 


in  (A. 3.1)  the  rotational  energy  can  be  written  as  a  scalar 


12  1  • (a)  -y 

t  I  w  =  -7  I  e  e ,  \ 

2  4  y  (a) 


1  *  (3.)  • 

Taking  the  above  for  the  Lagrangian,  L  =  ^  I  e^  e(a)  '  T^en 
9L  „  ....  ,_w,.  9L 


1  x  *V 

=  7  I  e 


.  .  r.  .  .  .  „yv  _  „  (a)[y  9~ 

and  defining  spin  by  S  =  2  e  ~ 7(a)" 


(A. 3. 3) 


gives 


.(a)  2  e(a)  . 

3ev  3e  V] 

Syv  =  I  V]  .  Since  .  =  0  ,  a  =  1 . 4  ,  the  components 

(a)  (a) 

=  -Sl;4  are  zero.  The  spatial  components  of  S*JJ  are  easily  seen 
from  (A. 3. 2)  to  be  just  the  components  of  the  angular  momentum,  e.g.. 


;12  „  I  !  e(a)l  42  _  1  (a)  2  ;1 

2  e  e(a)  2  e  6  (a) 


1  j  e  d£(i)2  ,IIe  dea)l 

2  e(i)l  dt  2  1  e  (i)  2  dt 


1  ^~fi) 

-  2  1  (~e(i)  X  "dt”^  =  1  W3 


■ 


.,'1 


.  .9  i  r  »  7  T - 


. 
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APPENDIX  FOUR 


Equation  (3.2.2). 


_  ra  pa 

0  “  saT|p  _  8ax,p  "  ap  8ax  xp  8aa 


0  (§aT  +  ^ax^lp  ^8ax  +  ^8ax\p 


-O1!  +  dr“  )(«  +  dg  )  -  <r“  +  ar“  > (g  +  d*) 


ap 


ap  °ax 


Jax 


xp 


xp  ua  a 


3a  a 


Subtraction  gives 


0  =  dg 


ax,p 


dra  g  -  ra  dg  -  dra  g  -  ra  dg 

ap8ax  ap  8ax  xp8aa  xp  8aa 


i.e. , 


o  =  (dg),  -  ar“  g 


3ax 


ap°ax 


df 


a 

xp8aa 


Noting  the  similarity  of  this  equation  to  the  first  equation  (which  implies 

r^p  =  —  (p  +p  -g  )),  cyclic  permutation  and  addition  will 

axg3p  2  k8ap,x  5xp,a  6ax,p"  y 

give 


1  ,,  P 

2  (dg  al 


+  dgh 


-  dg 

a  &ax 
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Equation  (3.2.3). 


pOt  _  pOt  pOt  ,  p  06  pT  p06  pT 

3y6  ~  36, y  3y,<5  xy  35  t6  3y 


In  geodesic  co-ordinates,  =  0  , 


_  jr®  j-pG 

dR  By«  '  drB«,Y  drBY,« 


2drB[fi,Yl  2drB[«lY]  ' 


,a 


drg6  is  a  tensor,  and  this  equation  holding  in  geodesic  co-ordinates 
implies  that  it  is  true  generally. 


Equation  (3.2.5). 


xg  ,r-P  „xg,,  P  1  ,  pN 

B  drK  =  B  (dg  ,  i  v  -  dg  ) 

p  ax  p  6  (a|x)  2  6ax 


...  N  1  -XG  ,  p  1  .TO  ,  p  .1  _XG  p 
( div)  —  —  B  i  dg  —  —  B  I  dg  +  —  B  dg 

k  2  p  x  6  g  2  p  g  6  x  2  p  ax 


(div)  +  i  C-BWP  -  B0TP  +  B0pX)  I  dgpa 


(A.  x  ,  1  /  T,T(ap)  (Gxp)  R(ap)T. 

(div)  +  ~2  (-B  -  B  +  B  )i  dgp(j 


(since  dgpQ.  is  symmetric).  Therefore, 


BTapdrPT  =  (div)  +  \  (-BX0P  -  BTpa  -  B0TP  -  BPT0  +  B0pX  +  BpOX)jx  dgpo 


=  (div)  +i  (Bp[OX]  +  B°[PX]  -  BX(pa)),  dg  „ 

2  x  pa 


- 


. 
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